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ABSTRACT 


The  present  report  develops  electrical  analogs  to  investigate 
multiple  cracks  on  simple  beams  and  more  complex  frame  structures. 

Analog  networks  provide  the  economic  tool  to  analyze  such  structures. 

The  effect  of  multiple  cracks  on  the  natural  frequencies  of  simple  struc¬ 
tures  is  studied  in  detail.  It  is  shown  that  closely  spaced  multiple 
cracks  are  indistinguishable  from  an  effective  single  crack.  A  severe 
crack  on  a  structure  can  be  identified  if  there  are  only  minor  cracks  in 
addition  to  the  major  one.  If,  on  the  other  hand,  there  is  more  than  one 
severe  crack,  then  the  damage  cannot,  in  general,  be  diagnosed  with  only 
three  frequencies  measurable.  Nevertheless,  a  minimum  number  of  cracks 
which  are  likely  to  be  present  in  the  structure  can  be  established. 

Characteristic  equations  are  developed  in  the  form  of  linear  systems 
for  cantilever  beam  and  general  frame  structures  with  multiple  cracks. 
Usefulness  of  relative-frequency-change  curves  are  demonstrated  and  rough 
guidelines  are  provided  to  aid  the  damage  diagnosis  process.  Several 
numerical  examples  are  included  which  illustrate  the  effect  of  multiple 
cracks  on  frequencies. 
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NOMENCLATURE 


Symbol 


e.j (*  L./L) 


Gi'  Gio 


*ii*  ^  i  2 


kc(kg) 


M1l’  Mi2 


P,  P’l 

Q,  Q*j 


crack  length 

half  depth  of  a  flexural  member 

location  of  the  crack  on  a  beam  with  only  one  crack 

distance  between  i-lst  and  ith  cracks  on  a  simple  beam; 
length  of  the  ith  beam  segment  in  a  frame  structure 

location  of  the  single  crack  which  is  equivalent  to  a  group 
of  cracks  (the  definition  of  equivalence  is  given  on  p.  24) 

modulus  of  elasticity 

voltage  sources  in  analog  T  circuit  for  the  ith  beam 
segment.  Equation  (10)  (no  subscript  i  for  a  one-member 
structure) 

conductances  in  analog  n  circuit  for  the  1th  beam  segment. 
Equation  (22)  (no  subscript  i  for  a  one-member  structure) 

L2/EIg2  (subscripted  when  referring  to  the  ith  segment) 

area  moment  of  inertia 

current  sources  In  analog  u  circuit  for  the  1th  beam 
segment,  Eq.  (23) 

total  number  of  cracks  on  a  structure 

total  number  of  column  (girder)  cracks  on  a  frame  structure 

length  of  a  beam;  characteristic  length  for  a  structure 

distance  between  i-lst  and  ith  cracks  on  a  simple  beam; 
length  of  the  1th  beam  segment  In  a  structure 

number  of  spans  of  a  frame  structure 

vector  of  unknown  moments 

resisting  bending  moment  at  the  i-lst  crack  in  a  cantilever 

resisting  moment  at  the  two  ends  of  the  ith  beam  segment 
(no  subscript  i  for  a  one-member  structure) 

defined  in  Equation  (66) 

variables  of  the  analog  n  circuit.  Equation  (20); 
subscripted  when  referring  to  the  1th  segment 


WWW 


Rj  relative  change  in  jth  modal  frequency  *  1  -  tOj/u»uj 

R.  relative  change  in  jth  frequency  due  to  an  equivalent  crack 

Jeq  causing  the  same  change  in  two  other  frequencies  as  the 
actual  damage 

S,$',T,T'  T  circuit  variables.  Equation  (7);  subscripted  for  the  ith 
segment 

U  matrix  of  resistances  or  conductances 

V.  ,V.  resisting  cross-shear  at  the  two  ends  of  the  ith  beam 

11  2  segment 

W,  W'  n  circuit  variables.  Equation  (20);  subscripted  for  the  ith 

segment 

X  coefficient  matrix  of  m  or  y' ,  Equations  (30),  (69) 

X<j>  submatrix  of  X  corresponding  to  the  jth  wall  of  the 
structure 

yi  transverse  deflection  at  the  ith  crack;  horizontal  dis¬ 
placement  of  the  ith  floor 

y.  ,  y-  transverse  deflection  at  the  two  ends  of  the  ith  beam 
1  2  segment  (no  subscript  i  for  a  one-member  structure) 

y  vector  of  transverse  deflections 

y^  rotation  of  the  ith  node  (prime  does  not  imply  derivative) 

rotation  of  the  node  r‘  Immediately  to  the  right  of  or 
yr  above  a  crack 

yii*  yi2  rotation  at  the  two  ends  of  the  ith  beam  segment 

Z-j>  zj0  resistances  in  T  circuit  for  the  ith  segment,  Eq.  (9) 

Z  coefficient  matrix  of  y,  Eqs.  (30),  (69) 

a  defined  in  Equation  (20) 

S  characteristic  value;  $k  =  w2pL4/(EI) 

6U  characteristic  value  of  the  undamaged  structure 

S( jth  characteristic  value 

jth  characteristic  value  for  the  structure  with  an  equiva- 
eq  lent  crack  having  two  other  characteristic  values  in  common 
with  the  actual  damaged  structure 


jth  characteristic  value  for  the  actual  damaged  structure 

characteristic  value  for  the  ith  beam  segment 
fracture  damage  =  a/b 

torsional  spring  constant  of  the  fracture  hing 
Poisson's  ratio 
lineal  mass  density 

sensitivity  number  for  the  ith  crack  on  a  simple  beam; 
sensitivity  number  for  the  crack  on  the  ith  member  of  a 
frame  structure 

sensitivity  number  for  a  crack  equivalent  to  a  group  of 
cracks 

jth  modal  frequency 

jth  modal  frequency  of  the  undamaged  structure 


1.0  INTRODUCTION 

The  study  of  damage  diagnosis  stems  from  the  need  to  assess  relia¬ 
bility  of  structures  which  have  been  subjected  to  unusual  levels  of 
excitation.  Such  a  diagnosis  problem  requires  the  availability  of  an 
analytical  model  and  of  measurements  of  some  structural  characteristics 
to  be  used  as  Inputs  to  the  model.  The  present  investigation  utilizes 
the  modal  frequencies  of  the  structure  in  detecting  and  identifying 
fracture  damage.  Damping  will  later  be  incorporated  into  the  model. 

Determination  of  the  crack  location  and  intensity  in  simple  struc¬ 
tures  with  only  a  single  crack  was  studied  in  [l]1.  There  It  was  shown 
that  a  crack  can  be  simulated  with  fracture  hinges  representing  the 
softening  effect.  The  classical  modal  shapes  for  the  beam  elements  wer 
assumed;  the  boundary  conditions  were  imposed;  and  the  resulting  set  o 
homogeneous  equations  were  solved  in  terms  of  the  modal  frequencies. 

Since  there  are  four  coefficients  associated  with  the  modal  shape  of  each 
beam  element,  the  order  of  the  resulting  system  is  four  times  the  number 
of  beam  elements  in  the  structure.  The  approach  soon  loses  its  feasibil¬ 
ity  as  the  structure  becomes  more  complicated.  The  present  report 
describes  another  approach,  namely,  the  use  of  electrical  analogy  in 
structural  analysis,  which  greatly  economizes  on  computation  of  modal 
frequencies  of  complicated  frame-structures. 

The  introductory  work  on  electrical  analogues  for  solving  static  and 
dynamic  problems  of  elastic  structures  is  surveyed  in  [2].  The  tech¬ 
niques  presented  there  are  suitable  for  analog  computer  applications.  In 
this  report,  a  modified  and  Improved  approach  to  free  vibration  is  devel¬ 
oped  which  is  readily  formalized  and  adaptable  to  digital  programming. 
Electrical  analogy  is  further  developed  to  simulate  cracks.  The  applica¬ 
tion  of  the  newly  developed  technique  is  extended  to  investigation  of  the 
behavior  of  multiple  cracks  on  structures.  When  the  method  is  applied  to 
free  vibration  of  multi-story  multi-span  frame  structures,  natural  freq¬ 
uencies  of  the  analog  circuit,  hence  modal  frequencies  of  the  frame,  can 

lumbers  in  brackets  denote  the  references. 
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be  determined  by  exciting  the  circuit  with  a  variable  frequency  excita¬ 
tion  and  by  establishing  the  excitation  frequencies  at  which  the  response 
is  maximum.  Analytically  this  corresponds  to  finding  zeroes  of  a  deter¬ 
minant.  The  order  of  the  system  for  a  multi -cell  structure  without  crack 
is  equal  to  the  number  of  nodes.  The  order  increases  by  two  for  each 
crack  present. 


2.0  BASIC  ANALOG  CIRCUITS 

The  mode  shapes  of  a  Bernoulli -Euler  beam  under  transverse  vibration 
are  given  by 

y(5)  3  A  coshBC  +  BsinhSC  +  C  cosBE  +  5  sins?  (I)2 


where  g4  *  po»2L4/£I;  p  is  the  lineal  density;  at  is  the  modal  frequency;  L 
is  the  beam  length;  (El)  is  the  beam  stiffness;  and  5  *  x/L  is  the 
normalized  axial  coordinate.  Four  variables  are  associated  with  each  end 
of  the  beam  element,  namely,  deflection,  slope  (or  angle  of  rotation), 
resisting  moment,  and  shear  force  (Fig.  I).  From  Equation  (1)  and  Its 
appropriate  derivatives  at  5  3  0,  the  variables  at  the  left  end  are 
obtained: 

yi  3  A  +  C 
y ! ’  3  (3  +  D)s/L 

Mj  3  (A  -  C)EIg2/L2  (2) 

V 1  3  (D  -  B)EIg3/L3 

Equations  (2)  can  be  solved  for  the  four  coefficients  In  terms  of  the 
four  variables 

A  3  (y!  +  hMx)/2 
B  3  (yi  -  hVj )L/2g 

C  3  (yt-  hMt)/2  (3) 

5  3  (yi  +  hV1)L/2g 


yz 
Yz 
M  2 

Vi  v2 

"inure  1.  Beam  Element  under  Free  Vibration, 
^lumbers  in  parentheses  denote  the  equations. 


where  h  3  L2  /Els2.  The  variables  at  the  right-end  of  the  beam  element 
are  obtained  from  Equations  (1)  and  (3)  with  5=1.  The  result  is 


y2  3  Ayx  +  ByiL/g  +  Cf^h  -  DVxhL/g  (4a) 
y£  =  Dy^/L  +  Ayl  +  BMxL/EIg  -  CVih  (4b) 
M2  =  Cyi/h  +  DyiEIg/L  +  AMX  -  BV^/s  (4c) 
V2  =  -By^/hL  -  Cy {/h  -  DM^/L  +  AVX  (4d) 

where 


A  =  (coshg  +  cosg)/2  B  =  (sinhg  +  sing)/2 

(5) 

C  =  (coshg  -  cosg)/2  D  ■  (sinhg  -  sing)/2 


Any  four  of  the  eight  variables  in  Equations  (4)  can  be  solved  for 
in  terms  of  the  remaining  four  variables.  The  mathematical  analogy 
between  electrical  circuits  and  beams  is  based  on  the  liner  transform 
(4).  The  type  and  properties  of  the  resulting  circuit  depend  on  the 
choice  of  the  independent  variables.  In  this  study,  angles  of  rotation 
are  analogous  to  voltages  and  resisting  moments  are  analogous  to  electri¬ 
cal  currents. 


2.1  T  Circuit 

When  it  is  imposed  that  sing  *  0,  slopes  and  shears  can  be  expressed 
in  terms  of  resisting  moments  and  deflections.  Thus,  from  Equations  (4) 


S  *  cothg  -  cotg 
T  *  cscg  -  cschg 


s'  *  cothg  +  cotg 
T'  *  cscg  +  cschg 


(7) 
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The  first  two  equations  In  (6)  can  be  rewritten  as 
yl  -  (Z  +  20 )MX  -  Z0M2  +  El 

(8) 

yi  »  Z0M!  -  (Z  +  Z0)M2  +  E2 

where 


Z  *  -  (S  +  T)L/2EIS 

Z0  -  TL/2EI0 

(9) 

Ei  =  (T ' y2  -  S * y ! ) 0/2L 

E2  *  (S'y2  -  T'y^s/ZL 

(10) 

Equations  (8)  are  the  Kirchhoff's  equations  for  the  active 
three  terminal  network  shown  in  Figure  2  (T  circuit)3  with  slope  and 
moment  (y * ,  M)  being  analogous  to  electrical  voltage  and  current, 
respectively.  The  quantities  Z  and  Ej  denote  resistance  and  voltage 
source,  respectively.  Negative  resistance  poses  no  difficulty  in 
analytical  and  numerical  analyses.  The  ends  of  a  beam  segment  are 
simulated  by  the  ports  of  the  circuit.  It  must  be  observed  that  the 
electromechanical  analogy  described  above  does  not  simulate  the 
differential  equation  of  motion,  but  instead  the  solution  based  on  the 
assumption  of  harmonic  motion  [2]. 

Boundary  conditions  at  the  ends  of  a  beam  element  can  be  simulated 
as  follows: 

a.  Free  end  or  simply-supported  end,  M  *  0:  the  corresponding  port 
of  the  circuit  to  be  left  open  {for  zero  current). 

b.  Fixed  end,  y'  *  0:  the  corresponding  port  to  be  short-circuited 
(for  zero  voltage). 

The  other  boundary  conditions  at  these  ends  and  the  last  two  equations  in 
(6)  are  utilized  to  express  the  voltage  sources.  Equations  (10),  in  terms 
of  the  moments,  as  will  be  illustrated  with  a  cantilever  beam  problem 
(Figure  3).  The  Kirchhoff's  voltage  law  applied  to  the  loop  yields 

(Z  +  Z0 )M !  +  Ex  *  0  (11) 


3Refer  to  the  appendix  for  the  development  of  the  circuit  equations. 
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where  Z  is  the  analog  resistance.  From  the  last  equation  in  (6)  with 

yi  «  M2  »  0, 

V2  *  0  -  (T * M 2  -  Sy2/h)8/2L  (12) 

from  which 

y2  «  HT'fVS  and  ll  -  T'2MlL/2EIflS  (13) 

The  second  equation  in  (13)  Is  obtained  from  Equation  (10). 

Substitution  of  Equations  (9)  and  (13)  Into  (11)  and  the  fact  that 
*  0  result  in 

T'2/S  -  S  *  0  (14) 

which  yields  the  characteristic  equation  of  a  cantilever  beam, 

1  +  coshg  cscg  »  0. 

In  derivation  of  Equations  (6),  It  is  theoretically  sound  to  Impose 
that  sins  *  0.  However,  during  the  numerical  search  of  the  natural 
frequencies,  sins  may  get  very  small  at  some  point,  or  the  structure  may 
actually  have  frequencies  at  or  near  sins  *  0.  We  shall  then  take  Into 
consideration  that  s  Is  directly  proportional  to  the  length  L  of  the  beam 
element.  Thereby,  this  problem  can  be  solved  by  dividing  the  beam  into 
two  sections  and  representing  each  section  by  an  Individual  circuit. 

This  will  be  Illustrated  with  a  simply-supported  beam.  The  beam  Is  arbi¬ 
trarily  divided  into  two  elements  with  a  length  ratio  of  4/6  (Fig.  4). 


Figure  4.  Simply-supported  Beam  Represented  with 
Two  Elements  and  its  Analog  Circuit. 
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With  reference  to  the  figure,  the  continuity  conditions  in  slope  and 
moment  (voltage  and  current)  at  the  common  boundary  is  preserved  under 
the  cascade  connection  of  the  two  basic  circuits.  The  loop  equation  then 
yields 

(ZL  +  Zjq  +  Z2  +  Z2(j)Mi2  ”  ^12  +  ^21  3  0  (15) 

Since  deflection  is  continuous  at  the  interface  (i.e.,  y12  0  y2i)»  can 
be  shown  that 


e12  3  Si *y i20i/2Li  E21  =  “  S2'yi262/21-2 

where  Li  *  .4L,  L2  *  .6L,  0i  *  .40,  and  02  3  *66  with 
From  Newton's  third  law  V12  *  V2i»  or 


pU2L4/EI. 


yi2  *  -  hM12(Si  +  Si)/(SX  +  S2)  (17) 

Equations  (9)  and  (15)  through  (17)  can  be  combined  to  yield  the  fre¬ 
quency  equation 

($! 1  +  S2')2  -  (Si  +  S2)2  *  0  (18) 

It  is  to  be  noted  that  for  the  particular  sectioning  in  this  illustra¬ 
tion,  the  same  problem  will  arise  at  the  fifth  natural  frequency  where  ^ 
=  2n  and  g2  3  3*.  Then  the  same  procedure  can  be  applied  by  further 
sectioning  the  beam.  Better  Initial  sectionings  are  of  course  possible 
than  the  one  presented. 

2.2  n  Circuit 

Similar  to  the  case  of  Equation  (6)  In  the  T  circuit.  Equation  (4) 
can  be  rewritten  such  that  resisting  moments  and  shears  are  expressed  in 
terms  of  deflections  and  slopes;  that  is, 

I  Mi  \  r  QL/S  WL/0  Q'  -W'  "I  / y 1 ' \ 

M  2  -WL/B  -QL/B  -W’  Q*  \y2‘ 

/  =  1  \  /  {19) 
VLl  Q‘  W'  P0/L  -P'B/L  /  y  1  V 

V2)  W'  Q'  P'b/L  -Pp/L  \y2  / 


P  =  (coshg  sing  +  sinhg  cosg)/<*  P'  *  (sinhg  +  s1ng)/a 

Q  =  (coshg  sing  -  sinhg  cosg)/a  Q*  ■  sinhg  sing/a 

W  =  (sinhg  -  sing)/o  W  =  (coshg  -  cosg)/a  (20) 

a  =  coshg  cosg  -1*0 


The  first  two  equations  in  (19)  can  be  rewritten  as 


Mj  a  (G  +  Go^l'  -  G0y2' 

+  Ii 

M2  *  G0y i '  “  (G  +  Go)y2' 

+  I2 

(21) 

where  G  and  I  are  the  analog 

conductance  and  current  source. 

respectively,  given  by 

G  =  (Q  +  W)EIg/L 

G0  =  -WEIg/L 

(22) 

II  ■  (Q'y 1  -  W'y2)/h 

I2  =  ( Q * y 2  *  W'y^/h 

(23) 

Equations  (21)  are  the  Kirchhoff's  current  equations  for  the  active  three 
terminal  network  shown  in  Figure  5  (n  circuit)4  with  G  and  I  denoting  a 
conductance  and  a  current  source,  respectively.  The  pair  (M,y')  are 
again  analogous  to  current  and  voltage.  It  must  be  noted  that  a  n 
circuit  can  be  derived  by  directly  inverting  the  system  in  (6).  However, 
the  resulting  derivation  requires  more  computational  effort  in  solution 
of  frame- structures.  As  an  illustration,  the  cantilever  beam  of  Figure  3 
will  now  be  simulated  with  a  n  circuit.  The  left  port  of  the  circuit  in 
Figure  5  is  shorted  yielding  the  circuit  of  Figure  6.  The  node  voltage 
equation  is 

(G  +  G0)y2'  “  I2  -  0  (24) 

The  last  equation  in  (19)  with  y!  =*  y L *  =  V2  *  0  yields 

y2  -  LQ,y2'/gP  (25) 

4Refer  to  the  appendix  for  the  development  of  the  circuit  equations. 


The  characteristic  equation  is  obtained  by  substituting  Equations  (22, 
23,  25)  into  (24).  Namely, 

Q  -  Q'2/P  =  0  (26) 


With  reference  to  (20),  Equation  (26)  Is  equivalent  to  1  +  coshecosg  *  0 
if  1  -  coshgcoss  *  0.  When  the  value  of  a.  Equation  (20),  approaches 
zero,  numerical  computation  diverges.  Similar  to  the  approach  in  T 
circuit,  the  beam  can  be  subdivided  into  elements  as  described  in  the 
previous  section. 

2.3  Simulation  of  Crack  with  Circuit  Analogy 

A  cracked  section  in  a  beam  is  modeled  following  the  method  of  frac 
ture  hinge  [1].  The  crack  is  mechanically  represented  by  a  torsional 
spring  of  spring  constant  *.  Slope  is  thus  discontinuous  at  the  cracked 
section,  the  discontinuity  being  given  by  &y‘  *  Mc/k.  Such  a  disconti¬ 
nuity  is  analogous  to  a  voltage  drop  In  the  circuit  theory.  The  crack 
can  then  be  simulated  with  a  resistor  of  resistance  -1/c.  Hence,  a  beam 
with  a  single  crack  Is  represented  by  two  circuits  (T  6r  n)  joined  by  a 
“crack  resistance,"  -1/tc.  Examples  will  be  given  in  the  next  chapter 
where  it  will  also  be  shown  that  crack  Intensity  is  quantified  by  the 
nondimensional  sensitivity  number  8  3  EI/ieL  [1]. 


3.0  MULTIPLE  CRACKS  ON  SIMPLE  BEAMS 


In  this  investigation,  damage  diagnosis  is  based  on  the  knowledge  of 
the  new  modal  frequencies  after  the  damage  has  occurred,  together  with 
the  knowledge  of  the  material  properties.  If  there  are  uncertainties 
about  the  material  properties,  then  the  frequencies  must  be  measured 
prior  to  the  damage  also.  This  point  will  be  explained  In  Section  3.2. 

Each  crack  in  a  structure  Is  characterized  by  two  variables  associ¬ 
ated  with  it,  namely,  location  (e)  and  Intensity  (y).  When  these  charac¬ 
teristic  variables  are  known,  changes  In  the  frequencies  can  be  computed 
with  the  method  presented  In  this  chapter.  This  will  be  called  the 
forward  problem.  The  inverse  problem  of  locating  and  identifying  the 
cracks,  however,  requires  (2k  +  1)  measurements  of  frequency,  where  k  is 
the  number  of  cracks  present.  The  extra  measurement  Is  needed  to  pin¬ 
point  the  cracks  because  of  the  multi ple-valuedness  of  the  equations. 

Only  in  some  special  cases,  2k  measurements  are  adequate  to  accomplish 
the  task.  On  the  other  hand,  the  number  k  Is  not  known  a  priori  in  a 
practical  situation.  Besides,  in  practice,  it  may  not  be  possible  to 
measure  more  than  the  first  few  modes.  In  application  It  is  not  possible 
to  locate  and  identify  the  cracks  deterministically  when  k  is  greater 
than  two.  However,  valuable  insight  may  be  gained  by  investigating  the 
forward  problem. 

This  chapter  studies  the  group  behavior  of  multiple  cracks  on  simple 
beams.  Guidelines  are  presented  for  damage  diagnosis  involving  multiple 
cracks.  Examples  are  given  to  illustrate  the  applications  of  these 
guidelines. 

A  beam  with  k  cracks  can  be  represented  by  k+1  circuits  joined  by 
resistors  simulating  the  cracks.  The  unknowns  in  the  T-circuit  analogy 
are  the  moments  (currents)  at  the  cracks  and  at  the  ends  of  the  beam. 
Hence,  there  are  (k+1)  unknowns.  On  the  other  hand,  the  unknowns  in  the 
H-circuit  analogy  are  the  slopes  at  the  ends  of  the  beam  and  at  the 
cracks.  There  are  two  unknown  slopes  at  each  crack.  Hence  the  order  of 
the  system  is  larger  with  n  circuits.  The  T-circuit  analog  is  found  to 


fm 


be  more  suitable  for  multiple-crack  analysis.  Therefore  this  section 
develops  the  characteristic  frequency  equation  for  the  general  case, 
establishes  the  conditions  under  which  multiple  cracks  become  equivalent 
to  a  single  crack,  and  Illustrates  the  Inverse  problem  when  k  *  1. 

3.1  Cantilever  Beam  with  Multiple  Cracks 

Figure  7  depicts  the  analog  circuit  for  a  cantilever  beam  with  k 
cracks.  Continuity  of  moments  at  the  cracked  sections  Is  preserved  via 
the  continuity  of  electrical  currents  through  the  "crack  resistors," 
-1/Kj.  The  order  of  the  system  of  mesh  current  equations  Is  (k  +  1). 

Namely, 

Ul  +  ZioJMi  -  Z10M2  +  En  =  0  \ 

“ZioMi  +  (Zj  +  Z10  +  Z2  +  Z20  ”  1/<i)M2  -  Z20M3  -  E12  +  E21  *  0  / 


~ZkoMk  +  (zk  +  Zko  +  zk+i  +  Zk+i,o  "  1/>ek)Mk+i  "  Ek2  +  Ek+i,o3  0 

From  Equations  (10),  the  analog  voltage  sources  {E^}  are: 

f  Eu  \  r T,-  0  ...  0  1  /y,  \ 

E12  Sj*  0  I  y2 

1  E21  I  ”S2‘  T2'  0  •  I  Y3 


_  B 

"  7U 


-T,'  S,’  0 


-k+1,1 


“Sk+i  Tk+i  \  yk+ij 


where  yj  Is  the  deflection  of  the  jth  cracked  section  and  yk+1  Is  the 
deflection  of  the  free  end.  These  deflections  are  solved  for  by  imposing 


Figure  7c.  T-circuit  Analog  for  the  Equivalent  Model  of  Figure  7b. 


the  shear  boundary  and  continuity  conditions.  That  is. 


V12  3  vi+i,i 

W  3  0 


i  =*  1,2,. ...,k 


where  the  first  subscripts  refer  to  the  numbers  of  the  elements  within 
the  beam  and  the  second  subscripts  1  and  2  denote  the  left  and  right  ends 
of  the  corresponding  elements.  Substituting  the  last  two  equations  in 
the  T-circuit  equations  (6)  into  (29)  and  rearranging,  we  obtain 

£  Z  y  =  X  m  (30) 

where  y  and  m  are  vectors  of  deflections  and  mesh  currents,  and  Z  and  X 
are  square  matrices  of  order  k+1  given  by 


Hence,  Z  is  a  tridiagonal  symmetric  matrix  and  X  is  an  upper  triangular 
band  matrix  with  a  band  width  of  3.  Solving  for  y  from  Equation  (30), 
provided  Z  is  non-singular,  we  have 

y  =  hZ-1X  m  (33) 

Substituting  (33)  into  (28),  and  then  the  result  into  (27)  along  with  the 
analog  resistances  from  Equations  (9) 

(X^X  -  U)m  =  0  (34) 

where 


U  *  0 


Sl+S2+20l6  T*2 


S2+S3+2020 


L°  ....  Tk  sic+sk+l+  2Vj 

with  a j  being  the  sensitivity  number  for  the  jth  crack  given  by  [1] 

0.  *  EI/Lk  .  =  3ir ( 1  -  v2)(r)/Y x[f(x)]2dx  (36) 

J  J  L  0 

where  y  =  a/b  is  the  relative  crack  depth  (Fig.  7a)  and  f(y)  is  the 
dimensionless  fracture  intensity  factor  for  symmetric  cracks  [1],  namely 

f (y )  -  (1  -  y)"1*5C1.122  -  2.363y  +  4.367y2  -  4.88y3  +  2.845y4  -  0.663y5] 

(37) 

As  Equation  (36)  implies,  e j  is  based  on  the  total  length  of  the  beam. 

For  non- trivial  solutions  of  (34),  the  determinant  of  its  symmetric 
coefficient  matrix  must  vanish,  yielding  the  characteristic  equation 


det  (XTZ-1X  -  U)  »  0  . 


When  the  crack  spacings  (e j )  and  sensitivities  ( e j )  are  known,  the 
natural  frequencies  can  be  computed  from  Equation  (38).  In  computing  the 
matrices  in  Equation  (38),  the  variables  S^,  T^,  Sj,  T}  (Equation  7)  per¬ 
taining  to  the  ith  beam  segment  are  computed  using  e-jg  =  (L-j/L) g 
where  8  is  based  on  the  total  length  of  the  beam. 

3.2  Damage  Diagnosis  with  a  Single  Crack 

Damage  diagnosis  using  the  present  model  is  accomplished  with  the 
knowledge  of  frequencies  after  the  damage  has  occurred.  Since  the 
characteristic  equations  of  structures  are  in  terms  of  dimensionless 
characteristic  values  {e> ,  the  material  properties  El  and  p  must  be  known 
accurately  to  compute  (8)  from  the  measured  frequencies  {<■>}.  If  the 
properties  are  known,  then  s  can  be  computed  from 


8  =  L(pu>2/EI)1/4 


(39) 


When  three  measurements  of  frequency  are  available  for  the  case 
k  =  1  (single  crack),  exact  location  (e)  and  severity  (e)  of  the  crack 
can  be  determined.  The  procedure  will  be  illustrated  with  the  cantilever 
beam  problem  developed  in  the  previous  section.  The  coefficient  matrix 
in  Equation  (34),  with  k  *  1,  is  of  the  form 


XTZ"1X  -  U 


I'll 

h12 


h12 

h22  "  208 


=  H(8,e,0) 


(40) 


where  the  8  values  in  Equation  (40)  are  the  post-damage  values  computed 
from  Equation  (39).  The  damage  characteristics  e  and  8  are  to  be  deter¬ 
mined  from  the  determinant  of  H,  which  can  be  written  as 


det[H(8,e,8)]  =  det[H(8,e,0)]  -  298hn  =  0  (41) 

where  e  *  l^/L  is  the  normalized  crack  location.  It  is  noted  that  when 
8  assumes  the  pre-damage  values  {gu},  det[H(s,e,0)]  is  equal  to  0;  that 
Is,  this  term  is  the  characteristic  equation  for  the  undamaged  (0  *  0) 
beam.  There  are  three  equations  emerging  from  Equation  (41)  for  the 
three  known  characteristic  values.  A  numerical  code  for  damage  diagnosis 


has  been  developed  in  which  e  is  varied  through  the  range  (0,1).  0  is 

computed,  for  a  given  e,  using  the  first  known  characteristic  value  in 
Equation  (41),  namely, 

9(e)  *  det[H(S(1),e,0)]/2B(1)hn  (42) 

where  the  superscript  on  g  denotes  which  characteristic  value  is  used. 

Then  with  the  second  known  characteristic  value,  g^, 

\ 

det[H(8(2),e,9)]  *  0  (43) 

where  the  value  of  0  is  now  substituted  from  Equation  (42).  A  zero* 
searching  routine  is  used  to  find  the  roots  of  Equation  (42)  which,  in 
general,  yields  multiple  solutions  for  e.  b^  can  then  be  used  to 
locate  the  crack.  With  e  known,  9  is  computed  from  (42). 

In  the  case  when  the  beam  deviates  somewhat  from  the  Bernoulli -Euler 

theory  or  when  the  material  properties  are  uncertain,  computing  {0}  from 

Equation  (39)  may  not  be  suitable.  The  knowledge  of  modal  frequencies 

prior  to  the  damage  is  then  necessary.  Since  <0  and  a>u  are  proportional 
2  2 

to  0  and  8^,  respectively,  where  u>u  and  bu  are  the  values  for  the  pre¬ 
damage  structure,  and  u  and  s  are  the  post-damage  values,  the  following 
relation  holds: 

B  =  8uU/a>u)1/2  (44) 

2 

It  is  assumed  in  (44)  that  the  proportionality  constant  between  uiu  and  b‘ 
does  not  change  after  the  damage  has  occurred.  The  characteristic  equa¬ 
tions  do  not  involve  any  material  properties  and  {0U}  can  be  computed 
from  the  characteristic  equation  for  the  undamaged  structure.  With 
(u»u)  and  {«))  measured  In  the  field,  {0}  can  then  be  computed  from 
Equation  (44)  and  the  above  procedure  can  again  be  employed  to  determine 
e  and  9.  Only  the  beam  length  L,  the  Poisson's  ratio  v,  and  the  slender¬ 
ness  ratio  (b/L)  are  required  to  determine  the  actual  crack  location  Lj 
and  the  crack  depth  a  (Eq.  36).  Since  v  is  nearly  the  same  for  most 


metals,  the  theory  becomes  Independent  of  the  specific  material  of  the 
structure  if  the  material  is  a  metal,  if  the  damage  is  presented  in  terms 
of  and  a,  and  if  {8}  are  computed  from  Equation  (44).  The  theory  is 
totally  independent  of  the  material  properties  and  dimensions  if  the 
damage  is  presented  in  terms  of  e  and  0.  The  data  are  more  conveniently 
presented  versus  the  relative  frequencies  U/a>u)  or  the  relative  fre¬ 
quency  changes  (1  -  u/a>u)  rather  than  versus  the  absolute  frequencies, 
whether  {8}  are  computed  from  Equation  (39)  or  (44).  In  this  investiga¬ 
tion,  relative  changes  are  chosen  as  the  means  to  convey  the  numerical 
results. 

In  practice,  the  response  of  structures  deviate  from  the  Bernoulli- 
Euler  theory.  It  is  therefore  more  accurate  to  compute  {8}  from  Equation 
(44).  Equation  (44),  in  effect,  adjusts  the  parameters  of  the  specific 
structure  such  that  the  structure  frequencies  match  those  predicted  by 
the  Bernoulli -Euler  theory. 

In  practice,  especially  for  structures  whose  frequencies  are  close 
to  each  other,  the  reduction  in  a  certain  frequency  due  to  damage  may  be 
small  while  the  reduction  in  the  next  frequency  may  be  large  enough  that 
the  post-damage  value  of  the  latter  frequency  may  drop  below  that  of  the 
former.  This  phenomenon  is  called  cross-over.  When  cross-over  occurs 
and  when  one  is  not  aware  that  it  has  occurred,  the  correspondence  estab¬ 
lished  between  the  pre-damage  and  the  post-damage  values  of  the  frequen¬ 
cies  will  be  in  error.  If,  in  such  a  case,  one  uses  Equation  (44),  which 
involves  the  ratio  of  the  pre-and  post-damage  values  of  the  frequencies, 
the  computed  characteristic  values  will  be  Incorrect.  On  the  other  hand, 
if  Equation  (39),  which  involves  the  structural  properties,  is  used,  a 
knowledge  of  the  pre-damage  frequencies  is  not  required.  Hence,  it 
becomes  Immaterial  which  pre-damage  frequency  a  certain  measured  post¬ 
damage  frequency  corresponds  to;  the  important  point  is  that  there  are 
some  frequency  values  available  which  satisfy  the  post-damage  character¬ 
istic  equation  (Eq.  41).  If  the  structural  properties  are  uncertain, 
they  can  be  Identified  by  measuring  one  frequency  prior  to  the  damage. 
That  is,  the  constant  K  in  the  relation  su  -  Ko^2  can  be  determined  by 
measuring  one  uu  and  computing  the  corresponding  8u  from  the  pre- 
damage  characteristic  equation.  The  same  constant  can  later  be  used  to 


compute  all  the  needed  post-damage  characteristic  values  b  from  the 
measured  post-damage  frequencies  at  by  means  of  s  ■  Kw1  .  In  this  case, 
however,  the  Rj  values  cannot  be  used  to  present  the  data,  since  the 
ordering  of  the  post-damage  frequencies  according  to  the  magnitudes  of 
their  values  does  not  correspond  to  the  pre-damage  ordering,  unless  the 
fact  that  cross-over  has  occurred  is  known. 

3.3  Equivalence  of  Multiple  Cracks  to  a  Single  Crack 

3.3.1  Approximate  Characteristic  Equations 

In  this  section,  we  study  the  conditions  under  which  the  identifica¬ 
tion  of  multiple  cracks  is  not  readily  distinguishable  from  that  of  a 
single  crack.  For  this  purpose, equivalence  of  two  cracks  on  a  simply- 
supported  beam  to  a  single  crack  is  first  established  analytically. 
Numerical  results  for  larger  numbers  of  cracks  are  then  presented.  It  is 
established  that,  when  equivalence  holds,  solution  to  the  inverse  problem 
of  damage  diagnosis  cannot  differentiate  between  single  and  multiple 
cracks. 

Figure  8a  depicts  a  simply-supported  beam  with  a  single  crack  and 
its  T-circuit  analog.  The  mesh  current  equation  together  with  the  Equa¬ 
tions  (9)  and  (10)  yield 

-(Sx  +  S2  +  29b)M  -  (S^  +  S2')y/h  -  0  (45) 

where  M  and  y  are  the  resisting  moment  and  deflection  at  the  cracked 
section,  respectively.  The  continuity  condition  V12  =  V21  at  the  crack 
allows  y  to  be  solved  for  in  terms  of  M.  The  characteristic  equation  for 
beam  with  a  single  crack  is  thus 

($!*  +  S2')2/(S1  +  S2)  -  (Si  +  S2  +  206)  =  0  (46) 


With  reference  to  Figure  8b,  the  characteristic  equation  for  the  case  of 
two  cracks  is  obtained  in  the  form  of  Equation  (38)  with 


"-Si'-Sj* 

v" 

Si+S2  T2 

= 

z  * 

.  V 

-S2*-S3,_ 

J  2 

(47) 


U  * 


Si+S2+29i8  T2 

T2  S2+S3+202b 


Now,  assuming  that  the  spacing  between  the~cracks  is  sufficiently  small, 
that  is,  e2  =  L2/L  =  e;  we  obtain  the  following  approximations: 


cosh  e2$  4  1  +  (es)2/2  ,  sinh  e2s  4  e$  +  (e8)/6 

cos  e2g  *  1  -  (es)2/2  ,  sin  e2fl  4  eB  -  (eb)3/6 

It  then  follows  from  Equations  (7)  that 

S2  4  2e0/3  ,  s2'  4  2/ eB  »  T2  »  eB/3  ,  T2'  4  2/eB 


(48) 


(49) 


Substituting  (49)  into  (47)  and  performing  the  necessary  operatidns,  we 
have,  for  the  matrix  of  coefficients  in  (34), 


XTZ_1X  -  U  4  X  j— ^  (St  +  S3) 

’  "  "  ’  d  (Us? 


■  i 

-i‘ 

+ 

hn 

h12 
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-i 

i 

h12 

h22 

1 

(50) 


where 


d  *  S1S3  +  2(S i  +  S3)(*b)/3  +  (*6)  /3 

hu  =  4(2  +  S ^S 3 ) / eB  +  CS[2S3  +  4Sj  -  S1S3(S1+  29l8) D 

h12  «  -2(4  +  sls3  +  S i S 3 ) / eB  “  2(S^  +  S3) 

h22  =  4(2  +  SiS3)/e6  +  CS i S3 2  +  4S3  -  S j S3 ( S3  +  29 2fi ) ] 


(51) 


The  above  analysis  Is  valid  at  points  where  Slt  Si',  S3,  S3'  are  of 
smaller  order  of  magnitude  than  (l/eB).  The  characteristic  equation  for 
the  two  crack  case  is  thus 
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(52) 


det(XTZ“1X  -  U)  *  j  (S/  +  S,’  )2 

“  -  -  "  UbT  ( 

-  (sx  +  s3 JCSi  +  S3  +  2(0!  +  02)0] 


+  O(l/E0)  =  0 


When  Equations  (46)  and  (52)  are  compared  and  It  is  noted  that  Slt  Sj', 
S2,S2'  in  Equation  (46)  are  nearly  equal  to  Sj,  Sx ' ,  S3,  S3',  respec¬ 
tively,  in  Equation  (52)  for  small  values  of  e2,  it  can  be  established 
that  the  difference  between  the  characteristic  equations  for  one-  and 
two-crack  cases  is  of  the  order  of  the  product  e2g;  that  Is,  the  differ¬ 
ence  is  larger  for  higher  modes.  The  effective  sensitivity  number  for 
the  two  cracks  approximates,  as  shown  in  (46)  and  (52),  the  sum  of  the 
individual  sensitivity  numbers,  0=0!+  02.  It  Is  shown  numerically 
that  closely  spaced  multiple  cracks  in  general  cannot  be  differentiated 
from  a  single  equivalent  crack. 


Similar  results  are  present  for  a  cl amped-c lamped  beam.  In  the  case 
of  a  cantilever  beam,  however,  the  expressions  are  more  complicated. 
Numerical  results  for  a  cantilever  beam  Indicate  the  existence  of  an 
equivalent  single  crack  to  closely  spaced  multiple  cracks.  Since  within 
10-15%  distance  of  the  free  end  of  a  cantilever  beam  frequency  changes 
are  small,  the  numerical  results  are  not  very  reliable  in  that  region. 

Equivalence  of  closely  spaced  cracks  to  a  single  crack  is  useful  in 
the  forward  problem  for  studying  the  behavior  of  multiple  cracks.  The 
implication  for  the  Inverse  problem  of  damage  diagnosis  Is  that,  using 
the  modal  method,  it  may  be  impossible  to  distinguish  between  closely 
spaced  multiple  cracks  and  a  single  crack.  In  special  cases,  the  equiva¬ 
lence  exists  for  some  configurations  of  multiple  cracks,  not  necessarily 
closely  spaced.  For  the  general  case,  we  may  establish  a  lower  limit  on 
the  spacing  of  cracks  which  is  assumed,  without  loss  of  generality,  to  be 
uniform. 


3.3.2  Lower  Limit  of  Spacing 

Figure  9  typically  Illustrates  the  lower  limit  of  crack  spacing  for 
a  cantilever  beam  when  uniformly  spaced  multiple  cracks  become,  as  a 


Dimensionless  Spacing  ( xl 0 3 ) 


Number  of  Cracks,  k 


Figure  9.  Lower  limit  of  crack  spacing  for  multiple 

cracks  on  a^ cantilever  to  be  distinguishable 
from  a  single  crack.  (All  cracks  are  equally 
spaced  in  each  case.  Agreement  in  the  third 
characteristic  values  is  =0.001.) 
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whole,  indistinguishable  from  a  single  crack.  In  Figure  9,  as  in  subse¬ 
quent  illustrations,  each  crack  assumes  a  sensitivity  number  e-f  of 
value  0.174  which  is  taken  to  be  the  lower  limit  for  failure  for  a  beam 
of  slenderness  ratio  of  0.05  (corresponding  to  a  relative  crack  depth 
y  =  a/b  of  0.6).  Figure  9  is  generated  by  the  trial  method.  With  the 
number  of  cracks  (k),  individual  sensitivities  (e-f )  and  the  first  crack 
location  (ex)  given  and  an  Initial  value  for  the  spacing  between  the 
cracks  assumed,  the  first  three  characteristic  values  s  are  computed  (the 
forward  problem).  These  values  of  8  are  called  the  actual  values.  The 
first  two  of  these  are  then  used  to  solve  for  a  single  crack  as  described 
in  section  3.2  (the  inverse  problem  with  k  3  1).  With  the  e  and  e  values 
for  the  equivalent  single  crack  thus  known,  the  third  characteristic 
value  for  the  equivalent  crack  is  computed  (the  forward  problem  with 
k  =  1).  The  8  values  for  the  equivalent  crack  are  called  the  equivalent 
characteristic  values,  the  first  two  of  which  are  the  same  as  the  actual 

(3)  (3) 

values.  If  the  difference  between  8gq  and  8^c^  is  not  equal  to  the 
error  allowance  (i.e.,  ±0.001),  then  a  new  value  for  the  crack  spacing  Is 
assumed  and  the  whole  procedure  is  repeated  until  the  required  agreement 
is  reached. 

There  are  three  cases  of  different  crack  locations  In  Figure  9.  In 

the  first  case,  the  first  crack  is  very  near  the  built-in  end,  at  a 

normalized  distance  of  ex  =  0.001  from  the  end.  In  the  other  two  cases, 

the  cracks  are  moved  toward  the  free  end  with  ex  *  0.1  and  ej  *  0.6, 

respectively.  In  all  the  cases,  (k  -  1)  equally  spaced  cracks  follow  the 

first  crack.  Table  1  shows  the  numerical  data  corresponding  to  the 

(3) 

ex  3  0.1  case.  If  the  tolerance  level  in  Sgq  Is  Increased,  the  curves 
in  Figure  9  will  shift  upward  (see  also  Figure  12).  As  the  number  of 
cracks  increases,  they  must  be  more  closely  spaced  to  be  representable  by 
a  single  crack.  The  lower  limit  of  spacing  depends  on  where  the  group  of 
cracks  is  located.  At  the  built-in  end  and  at  ex  3  0.6,  more  widely 
spaced  cracks  can  become  Indistinguishable  than  at  ex  3  0.1.  The  lower 
limit  of  spacing  is  dependent  also  on  e-j  values.  For  example,  when 
ex  3  0.1,  k  3  2  and  9x  3  02  a  0.01  (corresponding  to  y  3  0.2  for 
b/L  =  0.05),  the  smallest  crack  spacing  for  which  the  double-crack  damage 
becomes  indistinguishable  from  a  single-crack  damage  is  e2  *  0.025  under 


TABLE  1 

Numerical  Data  Corresponding  to  a  0.1  Case  In  Figure  9 


k 

Spacing,  e 

s131 

s(3) 

«l 

9eq 

®eq 

2 

.0065 

7.806 

7.807 

.3479 

.1032 

3 

.0033 

7.799 

7.800 

.5219 

.1032 

4 

.0021 

7.793 

7.794 

.6958 

.1031 

5 

.0015 

7.789 

7.790 

.8698 

.1029 

6 

.0011 

7.786 

7.787 

1.0438 

.1027 

|B<3)  -  £  0.001.  Figure  10  illustrates  the  relative  changes  in  the 

first  three  frequencies,  R.  =  1  -  o»./o»  •, corresponding  to  two  of  the 

J  J  **  J 

cases  In  Figure  9  where  aiuj  is  the  jth  undamaged  frequency.  Each 
change  Is  based  on  the  corresponding  ordinate  in  Figure  9.  For  any  given 
number  of  cracks  on  a  cantilever  beam,  the  largest  decrease  in  the  funda¬ 
mental  frequency  occurs  when  all  the  cracks  are  grouped  at  the  built-in 
end.  On  the  other  hand,  for  eeq  _>  0.3,  the  greatest  change  in  the 
second  frequency  Is  observed  at  a  location  of  0.55-0.60  (Fig.  11).  Such 
Information  can  be  utilized  to  set  lower  limits  on  the  possible  number  of 
cracks  on  the  structure  as  will  be  Illustrated  later. 

Table  1  lists  the  equivalent  locations  and  sensitivity  numbers  when 
cracks  are  grouped  at  e^  a  0.1.  The  table  indicates  that,  for  the 
equivalent  single  crack, 

9eq  -  0.174  k  ,  efiq  «  ex  +  e(k2~  1}-  (53) 

where  e  is  the  uniform  crack  spacing.  Numerical  results  indicate  that. 

In  general,  9eq  Is  approximately  equal  to  the  sum  of  individual  sensi¬ 
tivity  numbers  9i  and  that  the  equivalent  crack  location  eeq  is 
closer  to  the  cracks  with  larger  sensitivity  numbers. 
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Figure  11 


Relative  frequency  changes  for  a  cantilever  beam  with 
two  cracks  as  a  function  of  the  first  crack  location 
when  the  second  crack  is  kept  at  a  constant  distance 
from  the  first. 


3.3.3  Guidelines  for  Diagnosing  Multiple  Cracks 

In  general,  with  two  actual  characteristic  values  known,  more  than 
one  solution  exists  for  a  single  crack  as  was  pointed  out  In  Section 
3.2.  If  the  structure  does  indeed  have  only  one  crack,  the  crack  can  be 
diagnosed  accurately  with  the  knowledge  of  a  third  characteristic  value; 
that  is,  the  true  solution  can  be  pinpointed  among  the  set  of  discrete 
solutions.  In  other  words,  two  different  single-crack  configurations  can 
agree  in  at  most  two  characteristic  values.  On  the  other  hand,  a 
multiple-crack  configuration  can  have  the  same  three  characteristic 
values  as  a  single-crack  configuration.  Figure  12  compares  the  actual 
3(3)  values  due  to  two  cracks  on  a  cantilever  beam  with  b^3)  values 
which  would  stem  from  the  presence  of  either  of  the  equivalent  single 
cracks.  As  in  the  previous  illustration,  equivalence  is  defined  as 
8^  *  6 fo«"  j  *  1»2.  In  Figure  12,  the  first  crack  is  kept  at  a 
constant  location  (ex  a  0.1)  while  the  location  of  the  second  crack  is 
varied  relative  to  the  first.  In  the  figure,  e2  =  0  Implies  that  the  two 
cracks  coincide;  that  Is,  there  is  one  crack  at  e!  *  0.1  of  0  =  0.174 
+  0.174  *  0.348.  (0act*  J  *  1*2,3}  *s  computed  first  for  each  e2  value 
(the  forward  problem  with  k  *  2).  The  first  two  sact  values  are  then 

used  to  solve  for  the  equivalent  single  cracks  (i.e.,  to  solve  for  eea 

(3)  H 

and  eec|.  Bgq  is  then  computed  corresponding  to  each  (e»9)eq.  There  are 

two  single-crack  solutions  for  each  set  shown  by  the  and  0^ 

1  (3)  2 

curves  in  Figure  12.  The  second  solution  will  be  ignored  since  0  Is 

eq2 

(3) 

significantly  different  from  0^  over  most  of  the  range  of  e2.  For 

small  values  of  e2,  the  actual  and  equivalent  0^  values  agree 

closely.  In  such  a  case,  the  location  of  the  damage  can  be  established 

accurately  (eeq),  although  the  number  and  intensity  of  the  cracks  are 

uncertain.  (In  a  real  situation  only  {0^}  and  {0^,  efiC),  0^}  are 

known,  the  former  from  the  measured  {w-t}  and  the  latter  by  computa- 

J  (3) 

tion.)  As  the  spacing  between  the  cracks,  e9,  increases,  0V^  deviates 

eq  1 

more  from  the  actual  value  and  eeq  is  no  longer  indicative  of  the 
damage  location.  The  only  conclusion  that  can  be  drawn  in  this  case  Is 
that  there  is  more  than  one  major  crack  on  the  beam. 


[• 


[I 


Distance  between  the  Cracks,  ez 

Figure  12.  Third  characteristic  value  (S^)  variation  as  a  function  of  the 
location  of  the  second  crack  with  the  first  crack  fixed.  The 
other  curves  represent  (S)^  values  yielded  by  the  single  cracks 
-solved  using  the  first  two  actual  modes. 
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Variations  of  the  location  and  intensity  of  the  equivalent  crack 
(3) 

giving  rise  to  6^'  are  plotted  in  Figure  13.  For  small  values  of  e2, 
0eq1  *  0i  +  02  and  Seqi  “  ei  +  The  two  upper  curves  in  Figure 

12  intersect  a  second  time  at  e2  »  0.255.  This  means  the  following: 
two  cracks,  one  at  ei  =  0.1,  the  other  at  a  normalized  distance  of  0.255 
away  from  it  with  ex  =  e2  3  0.174  cause  the  same  decrease  in  the  first 
three  frequencies  as  a  single  crack  located  at  et  =  0.024  with  e  =  0.187 
(from  Fig.  13).  Furthermore,  there  may  be  other  combinations  of  multiple 
cracks  which  may  cause  the  same  frequency  changes.  It  is  clear  that 
probabilistic  methods  must  be  resorted  to  in  diagnosis  of  multiple  cracks 
when  an  insufficient  number  of  frequency  measurements  are  available. 
Nevertheless,  some  qualitative  conclusions  may  be  reached  with  the  help 
of  Figures  9-11  which  are  independent  of  material  properties.  In  partic¬ 
ular,  when  a  solution  has  been  obtained  for  a  single  crack  and  a  decision 
has  to  be  made  between  a  single  crack  a  group  of  closely-spaced  less 
severe  cracks,  one  may  argue  that  forma ciw..  of  a  single  crack  is  more 
likely  for  a  beam  undergoing  bending.  In  other  words,  a  small  mode  I 
type  crack  is  likely  to  propagate  under  bending  rather  than  other  cracks 
forming  nearby.  On  the  other  hand,  one  severe  and  several  minor  cracks 
may  exist  distributed  along  a  structure  in  which  case  the  objective  is  to 
be  able  to  diagnose  the  major  crack.  Guidelines  may  be  established 
qualitatively  based  on  the  frequency  change  curves  for  a  cantilever  beam. 
Figure  11,  in  which  case  we  may  propose  that 

a.  If  is  considerably  smaller  than  R2  and/or  R3,  then  the 
(major)  crack(s)  is  at  a  normalized  distance  greater  than  -0.45  from 
the  built-in  end  (Fig.  11).  In  addition. 


i.  if  R2  is  considerably  larger  than  R3,  the  crack  is  located 
around  0.45-0.65  relative  to  the  built-in  end  (Fig.  10b)(see  Example  3 
below) ; 

ii.  if  R3  is  larger  than  R2,  the  crack  is  at  a  distance  greater 
than  -0.7  or; 

iii.  if  R2  and  R3  are  comparable,  then  there  may  be  one  major 
crack  at  0.65-0.70,  or  two  major  cracks,  one  each  in  the  peak  regions  of 
the  R2  and  R3  curves. 
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b.  If  R2  Is  significantly  smaller  than  both  Rx  and  R3,  the  crack  is 
in  the  region  0.2-0. 3,  or  there  may  be  a  (major)  crack  at  0.2-0. 3  and 
another  (major)  one  at  a  distance  greater  than  0.8,  the  latter  being  more 
likely  if  R3  is  greater  than  Rx. 

c.  If  Rx  is  significantly  larger  than  both  R2  and  R3,  and  R2  and  R3 
are  rather  small,  the  crack  is  at  0. 1-0.2. 


d.  If  the  values  of  R^ ,  R2  and  R3  are  comparable,  several  possibil¬ 
ities  exist.  To  list  a  few,  there  may  be  a  number  of  cracks  distributed 
over  the  beam  (Examples  1  and  4  below);  the  crack  may  be  at  0.35-0.40;  or 
there  may  be  a  crack  at  0.0-0. 3  and  a  few  others  at  locations  greater 
than  -0.5. 

The  term  "major  crack"  employed  in  the  present  report  stems  from  the 
assumption  that  a  beam  fails  if  an  individual  crack  reaches  a  relative 
depth  of  0.6.  Whether  a  shallower  crack  may  be  considered  major  depends 
on  the  particular  application.  The  following  table  presents  eL/b  values 
corresponding  to  a  range  of  relative  crack  depths  y.  v  is  taken  to  be 
0.3.  Thus,  if  y  =  0.3  is  considered  serious  for  a  beam  of  b/L  a  0.05, 
then  a  crack  with  9  greater  than  0.024  would  signal  a  major  crack. 


9L/b 


.1 

.051 

.2 

.204 

.3 

.484 

.4 

.966 

.5 

1.825 

.6 

3.485 

In  the  present  study,  only  the  first  three  natural  frequencies  are 
assumed  measurable.  The  peak  locations  of  R2  and  R3  curves  in  Figure  11 
are  weak  functions  of  9  and  shift  to  the  right  as  9  incre** *..**.  For 
example,  as  9  is  increased  from  0.01  to  0.50  (effective  value  for  a  few 
closely  spaced  cracks),  the  peak  of  the  R2  curve  shifts  from  the  location 
0.530  to  0.563,  while  the  second  peak  of  the  R3  curve  shifts  from  0.710 
to  0.760.  (See  also  [1],  pp.  29-31.) 


Similar  guidelines  for  different  structures  can  be  established  based 
on  minimum  and  maximum  frequency  changes  although  the  procedure  may 
become  tedious  for  more  complicated  structures.  The  following  examples 
illustrate  some  of  the  applications.  In  all  the  examples,  a  slenderness 
ratio  of  b/L  s  0.05  is  assumed. 

Example  1.  To  illustrate  the  uncertainties  involved  in  damage  diag¬ 
nosis  when  the  three  Rj  values  are  comparable  in  value,  the  following 
relative  frequency  changes  are  assumed  to  have  been  computed  from  the 
measured  frequencies  of  a  cantilever  beam:  Rj  =  .0579,  R2  =  .0586, 

R3  =  .0591.  It  is  to  be  determined  whether  the  damage  mainly  consists  of 
one  major  crack. 

The  (actual)  characteristic  values  are  computed  from 


-  Vl/2 


where  {g^}  are  the  undamaged  characteristic  values.  The  first  three  of 

these  for  a  cantilever  are  1.8751,  4.6941,  and  7.8548.  From  Equation 

(3) 

(54)  in  particular,  a  7.6192.  The  first  two  actual  characteristic 
characteristic  values  yield,  according  to  the  procedure  given  in  Sec¬ 
tion  3.2,  one  solution  for  a  single  crack  located  at  e&q  =  .372  with 
9eq  =  .128  (corresponding  to  a  relative  crack  depth  of  y  =  .55).  These 
equivalent  values  are  then  used  in  the  forward  problem  with  k  3  1 

(Equations  31,  32,  35,  38)  to  determine  the  third  characteristic  value 

(3) 

due  to  the  equivalent  crack.  The  computation  thus  yields  g^  *  7.614. 
Then, 


(3)  .  fl (3)  _ 
act  0 


.0018  . 


(3) 

If  we  assume,  for  the  sake  of  argument,  that  g^j'  is  much  more  accurate 

sact  an<*  ^  tlie  accuracy  Bjct  1S  not  &etter  than  ±-0018,  the  equiva¬ 
lent  single  crack  solution  may  be  accepted  as  the  true  damage  diagnosis. 
On  the  other  hand,  the  given  Rj  data  were  actually  generated  in  a 


forward  problem  with  ten  cracks  of  equal  intensity  (8  was  taken  to  be  .01 
[i.e.,  y  *  .2]  for  all  the  cracks)  located  such  that  e1  =  .001,  e^  *  .094 
for  i  =  2, 3,..., 10.  Thus,  in  this  case,  even  though  the  actual  and 
equivalent  s^)  values  match  closely,  a  firm  diagnosis  cannot  be 
reached. 

Example  2.  The  frequency  changes  and  the  third  characteristic  value 

have  been  computed  from  the  measured  frequencies  similar  to  the  first 

(3) 

example:  Ri  *  .379,  R2  =  .210,  R3  =  .144,  sact  =  7.2672.  What  conclu¬ 
sions  on  the  damage  configuration  can  be  reached  from  these? 


The  first  two  characteristic  values,  computed  as  in  the  first' exam¬ 
ple,  yield  two  solutions  for  an  equivalent  single  crack.  One  of  the 

(31  (3) 

solutions  yields  a  'value  which  is  much  different  from  and  is 

(3) 

therefore  disregarded.  The  solution  which  gives  B  1  -  7.2678  Is 

e  =  .0077  and  9  .400.  The  largest  value  of  R,  for  any  number  of 

eq  eq  1 

9  =  .400.  The  largest  value  of  R*  for  any  number  of  cracks  occurs 

eq 

when  all  the  cracks  are  grouped  at  the  built-in  end  (Fig.  11).  From 
Fig.  10a,  when  there  are  two  closely-spaced  cracks  at  the  built-in  end 
with  81  3  92  ■  . 174,  Rx  is  equal  to  0.354.  For  three  cracks  in  the  same 
region,  R!  =  .432.  The  measured  value  of  Rx  given  above  is  .379.  Hence, 
it  can  be  concluded  that  there  are  at  least  three  cracks  on  the  beam 


(based  on  the  assumption  that  8  for  each  crack  is  not  greater  than  the 
limit  of  failure,  0.174,  for  b/L  *  .05).  In  principle,  a  large  number  of 
minor  cracks  (i.e.,  cracks  with  9  less  than  0.01  in  this  case)  can  cause 


the  same  reduction  in  the  first  three  frequencies  as  given  in  the  problem 
statement.  However,  it  may  be  argued  that  the  likelihood  of  having 


several  major  cracks  distributed  in  the  peak  regions  of  the  Rj  curves 
is  greater. 


Example  3.  This  example  illustrates  the  effect  of  minor  cracks  on 
Rj  values  when  minor  cracks  are  present  in  addition  to  a  major  crack  on 
a  cantilever  beam.  The  locations  and  intensities  of  the  cracks  are 
known.  The  characteristic  values  b^  are  determined  from  Equation  (38) 
in  Section  3.1  (the  forward  problem).  The  relative  frequency  changes  are 
then  determined  from  the  characteristic  values,  namely. 


m. 

F'S 


§. 


>  o> 
b\* 


V. 
.  ' 

r.V 

r.V 
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(55) 


a.  First  the  analysis  will  start  with  one  crack  (Fig.  14a),  i.e., 
k  *  1,  at  e  *  .54  (in  the  peak  region  of  the  R2  curve.  Figure  11)  with 
9  3  .174  (y  =  .6).  The  relative  frequency  changes  are  then  computed  as 
Ri  3  .029,  R2  =  .137,  R3  3  .014. 

b.  Now,  in  addition  to  the  crack  in  part  (a),  there  are  five  more 
cracks  of  intensity  9  3  .0026  (y  »  .1)  each,  distributed  on  the  cantile¬ 
ver  such  that  3  e2  3  e3  3  e5  3  .1,  e4  3  .24,  and  e6  3  .09  (Fig.  14b). 
The  actual  characteri Stic  values  are  computed  as  g  ^  3  1.840,  g^ 

3  4.351,  and  g(3)  =  7.768  from  which  Rj^  3  .037,  R2  3  .141,  R3  3  .022. 

If  these  values  of  Rj  were  computed  from  actual  measured  frequencies 
and  the  damage  were  to  be  diagnosed  (the  inverse  problem.',  then  (gU)} 
would  be  computed  from  Equation  (54)  in  Example  1.  When  this  is  the 

case,  and  g(2)  are  used  to  find  a  solution  for  a  single  crack 

which  is  (e__,  9__)  3  (.51,  .183).  This  single  crack  would  Droduce 
eq  cq 

,(3) 


3  7.845  or  R3eq  =  .0024. 
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Figure  14.  Cantilever  with  One  and  Six  Cracks 
(distances  are  normalized  based  on 
the  beam  dimension  in  the  corresponding 
direction) . 


c.  The  same  case  as  In  (b)  except  e-f  are  changed  to  .01  (y  =  .2) 
with  the  major  crack  remaining  the  same,  i.e.,  e4  ■  .174.  Hence,  the 
depths  of  the  minor  cracks  are  doubled.  The  values  of  e^  and  are 
used,  as  In  the  previous  cases.  In  the  forward  problem  to  compute 
Ri  *  .058,  R2  *  .151,  R3  *  .043.  Rj  and  R2  are  In  turn  used  In  the 
Inverse  problem  of  damage  diagnosis  with  k  *  1  to  obtain  (eeq,  0eq) 

*  (.47,  .224),  -  7.830,  R3eq  =  .0064. 

The  above  cases  indicate  that  the  value  of  R2  did  not  change  signif¬ 
icantly  when  minor  cracks  were  added  to  the  beam.  If  this  were  an  actual 
damage  diagnosis  problem,  In  which  case  only  {0^}  and  { R j}  would  be 
known  from  the  measured  data,  then  a  decision  would  have  to  be  made  about 
the  computed  (efiq,  8eq).  In  case  (b)»  R3  *  .022  whereas  R3eq  *  .0024. 

In  case  (c),  R3  *  .043,  but  R3gq  *  .0064.  Hence,  based  only  on  compari¬ 
son  of  R3gq  with  the  actual  values  of  R3,  the  solutions  for  (eeq,  9eq) 
would  be  rejected.  Nevertheless,  the  equivalent  damage  parameters  In 
both  cases  closely  identify  the  major  crack  located  at  the  distance  of 
.54  from  the  built-in  end  with  e  =  .174.  When  there  Is  no  minor  crack  on 
the  beam,  as  in  case  (a),  (eeq,  9eq)  would  be  computed  as  (.54,  .174). 

The  effect  of  the  minor  cracks  on  (e„,  eaJ  is  thus  seen  to  be  small. 
Hence,  It  is  concluded  that,  when  the  given  R.  values  exhibit  a  pattern 

J 

(3) 

such  as  in  this  example,  the  discrepancy  in  R3  and  R3eq  or  In  'ct 
(31 

and  can  be  Ignored  when  the  interest  is  in  diagnosing  the  major 
crack. 


fc: 


N* 


Example  4.  This  example  illustrates  different  crack  configurations 
which  cause  the  same  decreases  in  the  first  three  frequencies.  In  each 
case,  the  crack  locations  and  intensities  are  given  from  which  the 
characteristic  values  and  the  relative  frequency  changes  are  computed 
(the  forward  problem). 
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-r. 


The  three  crack  configurations  are 


a.  k  s  3,  e!  s  0.013,  e2  =  0.404,  e3  =  0.408 
01  s  02  =  83  3  0.05  (y  3  0.405); 

b.  k  =  4,  e3  =  0.001,  e2  =  0.15,  e3  =  0.23,  e4  =  0.25 
01  =  02  =  64  =  0.01  (y  =  0.2),  03  =  0.174  (y  =  0.6); 

c.  k  a  8,  ej  a  0.001,  e2  =  0.06,  e3  =  0.14,  e4  =  0.20, 
e5  =  e6  =  e7  =  0.10,  e8  =  0.19 

91  *  0.021  (y  3  0.28)  for  all  8  cracks. 

All  yield  the  same  changes  in  the  first  three  natural  frequencies, 
namely, 

Rx  *  0.099,  R2  3  0.108,  R3  3  0.089. 

There  may  likely  be  other  combinations  of  cracks  which  produce  the  same 
results.  When  the  first  two  characteristic  values  are  used  to  solve  for 
a  single  crack  (the  inverse  problem  of  damage  diagnosis),  ee£)  ■  0.384, 

0  3  0.249,  R3  3  0.081  are  obtained.  If  the  above  R-  values  had  been 

the  actual  measured  values  for  a  damaged  cantilever  and  if  the  true 
situation  had  been  as  in  (b),  then  the  major  crack  (03  =  .174),  which  is 
at  a  distance  of  0.381  from  the  built-in  end,  would  have  been  diagnosed 
accurately.  However,  uncertainties  always  exist  when  the  three  R-j 
values  are  comparable  to  each  other  as  illustrated  in  this  example. 

Example  5.  This  example  is  intended  to  show  that,  when  multiple 
cracks  exist  on  a  beam,  a  crack  which  is  located  in  the  peak  region  of  an 
R,  curve  affects  the  corresponding  frequency,  w .,  the  most. 

J  J 

a.  Given  k  *  1,  e3  3  0.001,  0!  3  0.174.  The  frequency  changes  are 
computed  as  (the  forward  problem) 

Rx  3  0.234,  R2  3  0.162,  R3  3  0.125. 

b.  A  second  crack  is  added  to  case  (a)  such  that  e2  3  0.546  (i.e., 
near  the  peak  region  of  the  R2  curve,  Figure  11)  and  e2  3  0.174.  Then 


c.  A  third  crack  is  added  to  case  (b)  such  that  e3  =  0.185  (l.e., 
near  a  peak  region  of  the  R3  curve)  and  q3  =  0.174  [parameter  values  in 
cases  (a-b)  are  preserved].  Then 

Ri  3  0.247  R2  *  0.301,  R3  =  0.250. 

If  03  3  0.01  instead  of  0.174,  then 

Rx  =  0.246,  R2  =  0.273,  R3  =  0.167. 

In  case  (b),  the  greatest  dcrease  relative  to  case  (a)  occurred  in  the 
second  frequency,  whereas  in  case  (c)  it  was  the  third  frequency  that 
suffered  the  largest  decrease  relative  to  case  (b)  although  this  decrease 
was  insignificant  for  q3  *  0.01.  The  example  illustrates  the  effect  of 
crack  location  and  intensity  on  the  frequencies  and  confirms  the  useful¬ 
ness  of  curves  such  as  Figure  11  to  diagnosis  of  damage  as  outlined  in 
the  guidelines  presented  earlier. 

Example  6.  In  all  the  above  examples,  0^  and  b^  have  been  used  in 

the  inverse  problem  of  diagnosing  the  major  crack  as  described  in  Section 

3.2.  It  has  been  found  numerically  that  using  the  pair  (s^,  (or 
12)  ( 3) 

( 6  ,  B  )  yields  better  results  in  some  cases  as  this  example 

illustrates. 

a.  Given  k  =  4,  ej  =  0.001,  e2  =  e3  =  0.15,  e4  ■  0.33, 

01  =  92  =  03  3  0.01,  04  *  0.174 

which  result  in 

R!  =  0.049,  R2  =  0.137,  R3  =  0.119 

b.  If  it  is  now  assumed  that  the  above  Rj  values  have  been 

obtained  from  the  measured  frequencies  and  If  {b^}  are  computed  from 
Equation  (54),  then  0^  and  can  be  used  in  the  inverse  problem  to 

solve  for  the  single  crack 


If,  instead,  s'  '  and  gv  ,  which  correspond  to  the  largest  two  of  the 
three  Rj  values,  are  used,  the  diagnosis  yields  two  possible  solutions, 
namely, 

eeq  '  °-372-  9eq  *  °*422’  Rieq  *  °*161 
and 

eeq  *  0.647,  0eq  =  0.217,  Rieq  =  0.015 

The  second  one  of  the  last  two  solutions  is  more  likely  since  Rleq*  0.015 
is  closer  to  the  measured  value  of  Riact  a  0.049.  Indeed,  the  location 
diagnosed,  egq  *  0.647,  is  close  to  the  location  of  the  major  crack  which 
is  at  the  distance  of  0.631  from  the  built-in  end. 

Other  examples  can  be  given  wherein  a  certain  pair  of  values 
yields  better  diagnosis.  However,  It  Is  not  clear  at  this  point  which 
pair  to  choose  in  each  case. 

It  should  be  noted  that  the  sensitivity  number  e  is  a  measure  of  how 
sensitive  the  natural  frequencies  are  to  given  crack  depth  and  location 
as  a  function  of  the  slenderness  ratio.  Decrease  In  natural  frequencies 
is  greater  for  larger  values  of  9.  Between  two  cracked  beams  with  the 
same  crack  location  and  relative  crack  depth,  9  for  the  more  slender  one 
will  be  smaller  giving  rise  to  smaller  Rj  values.  For  example,  for 
y  -  0.6  (and  v  *  0.3),  9  a  0.0697  when  b/L  =  0.02  and  9  =  0.0174  when 
b/L  =«  0.005.  Thus,  in  practice  it  is  relatively  harder  to  diagnose 
damage  in  slender  beams. 


ft 
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4.0  General  Two-Dimensional  Frame  Structure 

This  chapter  expounds  a  generic  model  of  a  two-dimensional  multi¬ 
storied  frame  (Fig.  15)  with  n  stories  and  m  spans  (thus  m  +  1  anchors). 

A  formalized  scheme  is  developed  for  obtaining  the  modal  frequencies  of 
such  a  structure  with  or  without  cracks.  In  the  T-circuit  analogy,  the 
characteristic  equation  for  the  structure  is  developed  through  the  mesh 
current  (moment)  equations.  Hence,  the  unknowns  are  the  moments  at  the 
frame  joints  and  at  the  cracks.  There  are  (p  -  1)  unknown  moments  at  a 
frame  node  where  p  beams  join  together.  Because  the  moment  is  continuous 
across  a  crack,  there  is  one  unknown  moment  for  each  crack.  It  follows 
that  the  order  of  the  resulting  linear  system,  which  is  equal  to  the 
number  of  unknowns,  is  [n(3m  +  1)  +  k]  for  a  frame  with  k  cracks.  On  the 
other  hand,  in  the  n-circuit  analogy,  node  voltage  (slope)  equations  lead 
to  the  characteristic  equation.  Slope  is  continuous  at  a  frame  joint; 
therefore,  there  is  one  unknown  slope  for  each  frame  joint  regardless  of 
the  number  of  beams  connected  there.  However,  there  are  two  unknown 
slopes  at  each  crack  location,  one  on  each  side  of  the  crack.  In  other 
words,  slope  is  discontinuous  across  a  crack.  The  structure  is  assumed 
rigidly  fixed  at  the  ground  level;  hence,  the  slopes  are  zero  at  the 
anchors.  The  order  of  the  system  is  thus  [n(m  +  1)  +  2k]  with  w  cir¬ 
cuits.  The  n-circuit  analogy  is  therefore  preferred  over  the  T-circuit 
analogy,  since  the  order  of  the  system  for  k  <  2mn  is  smaller  in  the 
former. 

The  procedure  to  obtain  the  characteristic  equation  of  the  frame 
structure  can  be  outlined  as  follows: 

1.  n  circuits  simulating  individual  columns  and  girders  joining  at 
right  angles  are  interconnected  such  that  the  boundary  conditions  (Ref. 

1,  p.  40)  at  the  analog-frame  joints  are  observed.  (It  will  be  shown 
that  a  network  diagram  need  not  actually  be  drawn.) 

2.  For  each  crack  present  on  any  element,  one  more  n  circuit 
(hence,  two  more  unknown  slopes-vol tages)  and  one  crack  resistance  (Sec¬ 
tion  2.3)  are  added  to  the  network.  A  girder  or  a  column  with  p  cracks 
is  thus  simulated  with  (p  +  1)  n  circuits  connected  in  cascade  via  the 
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crack  resistances.  The  first  and  the  (p  +  l)th  circuits  are  then  con¬ 
nected  with  the  rest  of  the  network  according  to  Step  1. 

3.  Kirchhoff's  current  law  is  written  in  terms  of  the  node  voltages 
which  correspond  to  the  angles  of  rotation  at  the  frame  joints  and  at  the 
cracks. 

4.  Under  free  vibration,  the  coefficient  matrix  obtained  in  Step  3 
must  be  singular.  If  the  damage  parameters  (i.e.,  location  and  severity) 
for  each  crack  are  known,  the  characteristic  values  of  the  structure, 

{ 0 } ,  can  be  determined  from  the  zero  determinant  of  the  coefficient 
matrix  (the  forward  problem).  If  the  characteristic  values  are  known  and 
a  diagnosis  of  the  damage  Is  required  (the  inverse  problem),  the  method 
in  Secton  3.2  is  used  to  find  a  solution  for  a  single  crack,  similar  to 
the  examples  in  the  previous  chapter. 

The  approach  taken  here  results  in  a  formal  scheme  which  can  be 
applied  to  frame  structures  without  referring  to  an  actual  network  dia¬ 
gram.  Reference  will  be  made  to  Figure  5.  The  word  "wall"  will  denote 
the  union  of  all  the  columns  on  the  same  vertical  line.  Each  node  and 
each  beam  element  is  identified  by  a  number.  Numbers  corresponding  to 
the  nodes  are  encircled  in  the  figure.  The  numbering  order  for  the  beam 
elements  begins  at  the  left  lowest  column,  proceeds  up  through  the 
columns  on  the  first  wall,  returns  back  to  the  second  anchor,  proceeds  up 
vertically,  and  continues  in  that  order.  Once  the  columns  are  finished, 
girders  on  each  floor  are  numbered  progressively  from  left  to  right 
starting  with  the  first  floor  and  continuing  on  with  the  upper  floors. 
Quantities  related  to  a  beam  element  such  as  G,  P,  Q,  etc.  are  sub¬ 
scripted  by  the  number  of  the  element.  Each  node  is  numbered  by  the 
column  under  the  node.  If  there  are  i  cracks  on  a  column  (girder),  the 
original  number  of  the  beam  element  refers  to  the  uppermost  (rightmost) 
segment  of  the  column  (girder)  which  is  now  represented  by  i  +  1  beam 
segments.  The  other  segments  are  numbered  following  the  largest  number 
in  the  scheme.  Due  to  discontinuity  of  slope  at  a  crack,  two  new  nodes 
are  created  on  the  two  sides  of  each  crack. 


4.1  Analog  Circuit  Equations  and  Boundary  Conditions 


Nodal  Equations: 

Figure  16a  shows  a  typical  node  r  at  the  intersection  of  column 
elements  r,  s  and  girder  elements  p,  q.  In  reference  to  Figures  5  and 
16a,  the  sum  of  the  branch  currents  entering  the  node  r  are  set  equal  to 
zero.  Namely, 

(yl  —  v 1 ) G  +  (y1.  -  y'  )G  +  (y,1  -  y')G  +  (y1  -  y')G 

vjr1  Jr'  po  3r'  qo  vjrk  Jr'  ro  yr'  so  (56) 

-y'(G  +G  +G  +G)+I„-I,+I--I,"0 
Jrv  p  q  r  s  p2  qi  r2  si 


where  (y^-  -  yp  is  the  analog  voltage  difference  between  the  nodes  i  and 

r  and  Gpo  is  the  conductance  which  connects  these  two  nodes.  The 
current  sources  appearing  in  Equation  56  are  dependent  on  the  transverse 
deflections  of  the  two  ends  of  the  beam  elements  to  which  they  corre¬ 
spond.  Under  small  deformation  theory,  vertical  deflection  of  the  nodes 
shown  in  Figure  15  is,  by  second  order  approximation,  equal  to  zero. 
Thus,  in  Figure  16a  yp2  *  yql  =  0,  yp2  and  yql  denoting  the  deflections 

of  the  right  end  of  the  pth  girder  and  left  end  of  the  qth  girder, 
respectively.  If,  in  addition,  the  nodes  i  and  j  do  not  bound  cracks, 
then  ypl  =  yq2  =  0.  Consequently,  Ipl  =  Ip2  *  Iql  =  Iq2  =  0  from 

Equation  (23).  On  the  other  hand,  Ir2,  for  instance,  is  not  zero  since, 
in  general,  yrl  *  0  and  yr2  =  y$1  *  0  (i.e.,  side  sway  is  allowed). 

Upon  substitution  of  Equation  (22)  into  (56), 

-^"p^p  *  VJq  +  "rA  *  "slV  '  yi"pl"p  ‘  yj"qlWq  '  yknrlMr 


-y;nsi“s  * 


(M  lI>2 ' 


V  +  rr2  '  's.>  ■  0 


where  npl  =  (El0/L)p^/(EIg/L)o  and  (EIs/L0)  is  the  characteristic  vali 


chosen  for  the  frame. 


junction  node 
(a) 


crack  nodes 


Figure  16 


(b) 


.  Typical  Nodes  on  the  Frame. 


In  the  formulation  of  crack  nodes.  Figure  16b  shows  two  nodes,  r  and 
r',  across  a  crack  on  a  beam  element.  Kirchhoff's  current  law  is  written 
for  the  node  r: 

(y‘  *  yr)Gro  +  lyr  “  yr){"K)  '  yrGr  +  lrt  s  0  {58) 

where  y '  is  the  rotation  of  the  node  r*  and  the  conductance  -ic  repre¬ 
sents  the  crack.  The  nodal  equation  for  the  node  r*  is  of  a  similar 
form.  Equation  (58)  is  valid  whether  the  crack  is  on  a  girder  or  a 
column.  Following  the  same  steps  which  led  to  Equation  (57)  we  may 
obtain 


■yr(nriGr  ”  of?)  ”  yr^8g^  ”  yqnri^r  +  ^ETgj0*r 


Assemblage  of  the  Nodal  Equations: 

From  the  small  deformation  theory,  points  on  girders  on  the  same 
floor  level  will  have  the  same  horizontal  displacement.  Namely, 


...  —  y. 


mn+i , j 


i  ®  l,...,n 


yij  3  yn+i , j  3  y2n+i , j  3  —  3  Wi.j  1  3  1 . n  *  j  3  1»2 

(60) 

Also,  yi2  =  yi+1  i  =  l,...,n-l.  On  the  other  hand,  deflection  conti¬ 
nuity  across  a  crack  implies  that  deflections  of  the  nodes  r  and  r'  in 
Figure  16b  are  equal.  Thus,  there  are  (n  +  k)  unknown  deflections 
associated  with  the  frame,  k  denoting  the  total  number  of  cracks  on  the 

frame.  Designating  by  y. (y^  =  y^2,  i  =  l,...,n)  the  horizontal  displace¬ 
ment  of  the  ith  floor,  we  obtain  an  (n  +  k)  vector  of  displacements  and 
an  [n(m  +  1)  +  2k]  -vector  of  rotations  as 

y  5  {y !  y2  •••  yn  yn(2m+l)+l  •••  yn(2m+l)+k^T 

l'  S  ***yn(m+l)  yn(2m+l)+l  yn(2m+l)+l  *  * *yn(2m+l)+k^T  (61) 


where  y n ( 2m+l ) -*-1  des*9nates  the  transverse  deflection  at  the  ith  crack, 
and  yn(2m+i)+i  and  yn(2m+l)+i  des-i9nate  the  rotations  immediately  to  the 


left  of  (below)  and  to  the  right  of  (above)  the  ith  crack  which  Is  on  a 
girder  (column).  The  unknown  deflections  can  be  related  to  the  rotations 
via  the  kinetic  equations,  defining  the  sideway  motion  of  the  floors,  and 
the  shear  continuity  conditions  at  the  cracks.  Cross  shears  at  the 
column  ends  are  the  axial  driving  forces  for  the  floors  [1,  p.  42].  The 
axial  acceleration  of  a  floor,  which  is  composed  of  all  the  girders  on 
the  same  level,  is  the  same  as  the  transverse  acceleration  of  the  column 
ends  at  the  nodes  where  the  columns  join  the  floor.  Thus  for  the  ith 
floor  (Figure  (17) 


m 


Erv 

j=ol 


p+m-1  q+k-j-l 

[Vjn+i,2  Vjn+i+l,l]  ~  j 

2  PjLj  +  2  pjLj 

j=p  J  J  j=»q  J  J 

(62) 


where  p  =»  (m  +  l)n  +  (i  -  l)m  +  1  in  accordance  with  the  numbering  order 
described  earlier  and  ki  is  the  total  number  of  cracks  on  the  ith 
floor.  There  is  one  such  equation  for  each  floor.  The  term  in  the 
braces  on  the  right-hand  side  of  Equation  (62)  is  the  total  mass  of  the 
ith  floor.  Equation  (62),  upon  substitution  of  the  last  two  of  Equations 
(19)  for  shears,  will  relate  deflections  to  rotations.  Figure  18b  shows 
the  diagram  of  a  frame  with  m  »  1,  n  ■  2,  k  *  2.  There  are  eight  nodes, 
eight  beam  elements  (hence  eight  interconnected  n  circuits)  and  four 
unknown  deflections.  The  first  floor  has  one  crack  on  it.  The  kinetic 
equation  for  the  first  floor  with  pg  =  P5,  has  the  form 


V12  +  V32  “  V21  ’  V71  =  «2y  105^5  +  Lg) 
where,  for  instance. 


(63) 


V32  =  ^^31^3  “  P 3y  lS3/L3  )  (64) 

V71  3  -^CQ7y3  +  W7y7  +  (P7yj  -  P7y7)s7/l-7] 
from  Equation  (19)  with  y31  =  y31  *  0,  yx  =  y32  =  y7i, 

I  1  1  II 

y3  =  y32  38  y7i,  y?2  s  Y7t  y?z  s  Y7*  Substituting  shear  expressions. 
Equations  (64),  and  the  similar  equations  for  V12  and  V2l,  into  Equation 
(63)  and  rearranging,  we  obtain 
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n+i  +1 


m 


^  *L‘^0{  13f>  1  +  n23P2  +  n33p3  +  n73P7  +  ^53(85  +  68^1 

-  023 P “  ^73? 7y7^  =  ( n 12^ 1  ■  n22Q2^yi  "  n22^2y2 
+  (*>32^3  “  072^7^3  "  n72^7y7  (65) 

where 

nu  3  (EI)t  (S^/^  )J/C(EI)  0(S0/L0)J3  •  (66) 

The  kinetic  equation  for  the  second  floor  can  be  developed  similarly. 
Application  of  shear  continuity  at  a  crack  will  now  be  demonstrated  for 
the  crack  on  the  first  floor  {Fig.  18b).  Hence  V82  3  V5l,  or  from  the 
last  two  of  Equations  (19) 

4  (w8yi'  +  Q8y8’  "  Pay^s/Ls)  3  4  (Qsy8  +  W5y3  +  PsysBs/Ls)  (67) 

M8  "5 

since  y81  3  y52  3  0.  The  following  expression  is  obtained  by  rearranging 
Equation  (67)  and  noting  that  h5  *  h8  and  B5/L5  *  b8/L8: 

(fj)(ps  +  p8)y8 s  way  1  -  wjyj  +  Qays  ■  Qsys  (68) 

For  the  general  frame  structure,  there  are  n  kinetic  equations  of 
the  form  of  Equation  (65)  and  k  shear-continulity  equations  of  the  form 
of  Equation  (68).  These  (n+k)  equations  can  be  arranged  in  a  matrix  form 
to  solve  for  the  deflections  in  terms  of  the  rotations,  namely, 

y  3  (t)  z~1x  y'  •  (69) 

\p/o  -  -  - 

On  the  other  hand,  the  nodal  equations  (57)  and  (59)  can  be  arranged  in 
the  form 


where  I  Is  the  vector  of  current  sources  and  is  related  to  the  deflec¬ 
tions  via  Equations  (23);  that  is, 

•  m) 

After  Equation  (69)  is  substituted  into  Equation  (71),  then  Equation  (71) 
into  Equation  (70),  and  it  is  noted  that  HjHj  =  XT,  the  slope  equation 
results; 

(XTZ-1X  -  U)y *  =  0  .  (72) 


This  is  the  general  form  of  the  result  for  any  planar  frame  with  any 
number  of  cracks  on  it.  The  coefficient  matrix  in  (72)  Is  symmetric. 
The  individual  matrices,  in  general,  have  the  forms 
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(73) 


(74) 


with  y  and  y'  arranged  as  in  Equation  (61) 


where  the  subscripts  c  and  g  refer  to  cracks  on  the  columns  and  on  the 

girders,  respectively,  with  k  and  k  being  the  total  number  of  cracks  on 

c  g 

the  columns  and  on  the  girders  (k  *  k  +  k).  It  should  be  noted  that 
(i )  c  9 

X  ,  U0,  and  Z0  for  the  case  of  no  crack  are  modified  when  cracks  are 

introduced  on  the  frame.  The  matrices  U  and  Z  are  symmetric.  A  formal 

procedure  will  now  be  described  to  establish  the  above  matrices  without 

the  need  to  draw  the  actual  analog  circuit. 

4.2  Procedure  to  Establish  the  Matrices 

In  the  following,  for  simplicity  of  explanation,  properties  of  the 
beam  elements  are  assumed  uniform  throughout  the  frame, -that  is,  n^  *  1 
for  all  i,j,  except  for  the  lengths.  At  the  end,  corrections  for  the 
general  case  are  mentioned. 

1.  To  establish  the  U  matrix: 

U  is  composed  of  the  coefficients  of  the  rotations  y^ '  in  equa¬ 
tions  such  as  (57)  and  (59).  The  ith  row  In  the  matrix  stems  from 
Kirchhoff's  current  law  written  for  the  ith  node,  the  first  n(m+l)  rows 
being  for  the  nodes  on  the  frame  joints  and  the  next  2k  rows  for  the 
nodes  at  the  cracks.  Hence, 

a.  the  diagonal  entry  u^  will  be  the  sum  of  the  Q  values. 
Equation  (20),  of  the  beam  elements  adjoining  at  the  node  for  which 
Kirchhoff's  current  law  is  being  written.  If  this  node  is  at  a  crack, 
then  u-j-j  is  given  by  the  Q  value  of  the  element  ending  at  this  node 
minus  the  M  value  for  the  corresponding  crack,  where  M  *  1/se.  Diagonal 
entries  of  Bc  and  will  thus  be  of  the  form  (Qt  -  Mr).  And, 

b.  the  off-diagonal  entry  u^.  will  be  W^,  Equation  (20),  if  the 
nodes  corresponding  to  the  columns  i  and  j  of  the  matrix  are  linked 
directly  by  the  kth  beam  element.  If  these  two  nodes  are  linked  by  a 
crack,  i.e.,  they  are  the  nodes  on  the  two  sides  of  a  crack,  then  u^j 

is  given  by  the  M  value  of  the  crack.  Otherwise  u.  .  ■  0.  Due  to 

•  J 

symmetry,  u^.  =  u ^ .  The  ith  row  of  U0  In  Equation  (74)  is  unchanged 
relative  to  U0  of  the  no-crack  (k=0)  case  if  the  1th  node  is  not  adjacent 
to  a  cracked  element. 
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2.  To  establish  the  Z  matrix: 

The  matrices  Z  and  X  stem  from  the  kinetic  equations  and  shear  condi¬ 
tions.  The  first  n  rows  in  them  represent  the  kinetic  equations  govern¬ 
ing  side  sway  of  the  n  floors.  The  last  k  rows  correspond  to  shear 
continuity  conditions  at  the  cracks.  Z  consists  of  the  coefficients  of 
deflections  y,  and  its  entries  can  be  generated  as  follows: 

a.  The  diagonal  entry  z^  for  i=l,2,...n  will  be  the  sum  of 
the  0^  values  of  all  the  girder  elements  on  the  ith  floor  plus  the  sum 
of  the  P  values.  Equation  (20),  of  all  the  columns  adjoining  at  the  1th 
floor.  If  a  column  adjacent  to  the  ith  floor  is  cracked,  then  only  the  P 
value  of  the  column  segment  nearest  to  that  floor  will  be  included.  0^ 
can  be  expressed  in  terms  of  the  characteristic  value  0O  as  0k  =  ek0o. 

If  the  frame  is  homogeneous  in  properties,  then  ek  =  Lk/L0.  Diagonal 
entries  z,,  for  the  last  k  rows  (i.e.,  diagonal  entries  of  the  D  matri- 

11  m 

ces)  will  be  the  sums  of  the  P  values  of  the  two  beam  segments  on  the  two 
sides  of  the  corresponding  crack. 

b.  The  off-diagonal  entry  z^  i+1  of  Z0  will  be  the  negative  of 
the  sum  of  the  P'  values.  Equation  (20),  of  all  uncracked  columns  joining 
ith  and  (i+l)th  floors,  z^  ^  *  0  if  all  the  columns  between  those 
floors  are  cracked,  z^j  *  0  for  n  £  j  £  i+2  and  z^ ^  *  z^.  Z0  is, 
hence,  a  tri diagonal  symmetric  matrix.  If  any  one  of  the  columns  adja¬ 
cent  to  the  ith  floor  is  cracked,  then  z^.  for  1  <  i  <  n,  n  <  j  £  n  +  kc 
(i.e.,  entries  of  C_),  will  be  zero  except  for  that  (those)  j  value(s) 
which  correspond(s)  to  the  deflection(s)  at  the  crack(s)  nearest  to  the 
ith  floor,  in  which  case  z^j  is  equal  to  the  negative  of  the  P'  value 
for  the  column  segment  which  links  the  ith  floor  to  that  crack.  If  none 
of  the  columns  adjacent  to  1th  floor  is  cracked,  then  ith  row  of  C  is 

**  w 

zero  and  the  ith  row  of  Zq  is  unchanged  relative  to  Zq  of  the  no-crack 
(k=0)  case. 

zi  i+1  for  n  <  1  1  +  kc  (entries  of  Dc)  will  be  zero  if  there  is  no 
other  crack  between  the  (i-n)th  column  crack  and  the  floor  level  above 


this  crack.  Otherwise,  ^  will  be  equal  to  the  negative  of  the  P* 
value  for  the  column  segment  which  links  the  (i-n)th  and  the  (i-n+l)st 
cracks. 

z.  for  n+k„  <  i  <  n  +  k  (entries  of  D  )  will  be  zero  If  there 
i,i  +1  c  —  g  »g 

is  no  other  crack  between  (i-n-kc)th  girder  crack  and  the  wall  to  the 
right  of  this  crack.  Otherwise  zjj+i  will  be  equal  to  the  negative 
of  the  P'  value  for  the  girder  segment  which  links  the  (i-n-kc)th  crack 
and  the  crack  on  its  right. 


Hence,  Dc  and  are  diagonal  if  at  most  one  crack  exists  on  each 
beam  element.  If  more  than  one  crack  exists  on  any  column  or  girder, 
then  D£  and/or  Dg  are  accordingly  tridiagonal  symmetric  matrices. 

3.  To  establish  the  X  matrix: 


X  consists  of  the  coefficients  of  rotations  y . 1  in  equations  such  as 

(65)  and  (68).  A  square  submatrix  x(k)  in  Equation  (73)  corresponds  to 

rotations  of  the  frame  joints  on  the  kth  wall  (Fig.  15). 

(k) 

a.  The  diagonal  entry  x^  ,  which  stems  from  the  kinetic  equa¬ 
tion  for  the  ith  floor,  will  be  given  by  the  Q'  value.  Equation  (20),  of 
the  column  (column  segment)  under  the  frame- joint  node  n(k-l)+1  minus  the 


Q  value  of  the  column  (column  segment)  above  the  same  node 

( k) 

no  cracked  column  on  the  kth  wall,  then  x:*  =  Q',.,  ....  - 

/ 1.  \  11  iHK-l/'rl 

for  i=l,...,n-l  and  x^;  =  Q'nk  . 


.  If  there  is 
Qn(k-l)+i+l 


iv.) 

b.  The  off-diagonal  entry  x}  ^+1  =  "Wn(k-l)+i+l  ^that  1S»  the 
negative  of  the  W  value.  Equation  (20),  for  the  column  above  node 
n(k-l)+i)  for  i=l,...n-l  if  the  column  of  the  kth  wall  between  the  ith 


and  (i+l)th  floors  (i.e.,  the  column  above  node  n(k-l)+i)  has  no  crack  on 
(k) 

it.  Otherwise  xj  i+1  =  0,  and  that  entry  on  the  ith  row  of  Fc  which 
corresponds  to  the  rotation  of  the  crack  node  nearest  to  and  above  the 
node  n(k-l)+i  will  be  equal  to  the  negative  of  the  W  value  of  the  column 

segment  linking  the  node  n(k-l)+i  and  the  said  crack  node.  x)V  =  0  for 

(k)  (k)  (k) 

j  _>  i+2  and  xiV  =  -  xv.V  for  i  *  j.  (Entries  of  each  V  are  numbered 

I  J  J  I  • 

(k) 

independently  of  the  other  submatrices.)  X  is  thus  diagonal  if  each 


column  of  the  kth  wall  has  at  least  one  crack.  Otherwise  it  is 


tridiagonal.  is  unchanged  relative  to  of  the  no-crack  (k=0) 

•m 

case  If  there  is  no  crack  on  the  kth  wall. 

If  the  column  below  node  n(k-l)+i  is  cracked,  then  that  entry  on  the 
ith  row  of  Fc  which  corresponds  to  the  rotation  of  the  crack  node  nearest 
to  and  below  the  node  n(k-l)+i  will  be  equal  to  the  W*  value  of  the 
column  segment  linking  the  two  nodes.  Except  for  this  and  the  above 
mentioned  cases,  entries  of  Fc  are  zero. 

c.  The  E  matrices  in  Equation  (73)  couple  the  crack  nodes  with  the 
frame- joint  nodes.  Rows  n+1  through  n+k  of  X  (i.e.,  E  and  H  matrices) 
are  filled  in  as  follows:  The  entry  on  the  matrix  column  which  corre¬ 
sponds  to  rotation  of  the  node  on  the  left  (or  lower)  side  of  the  crack, 
at  which  the  shear  continuity  condition  is  being  written,  is  equal  to  the 
Q‘  value  of  the  girder  (or  column)  segment  which  links  this  node  to  the 
one  on  its  left  (or  below  it).  The  entry  corresponding  to  the  node  on 
the  right  (or  upper)  side  of  the  crack  is  equal  to  minus  the  Q'  value  of 
the  segment  linking  this  node  to  the  node  on  its  right  (or  above  It). 
These  two  entries  are  within  the  subroatrlx  Hn  (or  H  ).  The  entry,  corre- 
spondlng  to  the  node  on  the  left  of  (or  below)  the  left  (or  lower) 
crack-node,  is  equal  to  the  W‘  value  of  the  segment  linking  these  two 
nodes.  Finally,  the  entry  corresponding  to  the  node  on  the  right  of  (or 
above)  the  right-side  (or  upper)  crack-node  Is  equal  to  the  negative  of 
the  W*  value  of  the  segment  linking  the  two  nodes.  These  two  entries  can 
be  within  the  submatrices  E  or  H  depending  on  whether  there  are  one  or 
more  cracks  on  a  girder  (column).  Ec  k  *  0  if  there  is  no  crack  on  the 
kth  wall;  En  »  0  if  there  Is  no  crack  on  any  of  the  girders  adjacent  to 
the  kth  wall . 

4.  Corrections  for  the  general  case: 

If  beam  properties  are  not  uniform  throughout  the  frame  (each  column 
or  girder  still  has  uniform  properties  within  itself),  then  the  following 
corrections  are  necessary: 

a.  Multiply  each  term  In  the  U  matrix,  except  the  terms  of  the 
form  1/80^,  by  the  corresponding  n!  value,  that  is,  Qp  and  Wp  by  npl. 
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b.  Multiply  each  term  in  the  first  n  rows  of  the  X  matrix  by 

•m 

the  corresponding  n2  value,  that  is,  Qp‘  and  W'p  by  np2. 

c.  Multiply  each  term  in  the  first  n  rows  of  the  Z  matrix  by 
the  corresponding  n3  value,  that  is,  Pp  and  Pp'  by  np3  and  er  by  nr3. 

d.  Multiply  the  ith  row  of  Z,  n+1  £  i  <_  n+k,  by  ( 0^/L^)/ ( ) 
where  t  is  the  number  of  the  element  on  which  the  crack  for  which  the  ith 
row  represents  the  shear  continuity  condition  is  located. 

The  above  procedure  will  now  be  illustrated  with  an  example. 

Example  1.  The  illustrated  frame  is  depicted  in  Figure  18. 
a.  Frame  without  cracks  (Fig.  18a)  with  different  beam  properties. 
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where,  for  convenience,  the  rotations  and  the  deflections  at  the  nodes 
are  written  above  the  corresponding  columns  of  the  matrices.  The 
det  (XTZ*1X-U)  then  yields  the  modal  frequencies  of  the  frame  without 
crack. 


b.  Frame  with  two  cracks 

(Fig.  18b).  The 

frame  is  now  assumed 

uniform.  Hence,  n^ 

1 
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(79) 

are 

the  same  as  those  in  the  second  row  in  Equation  (76)  (except  for  the 
factors  npl  .  This  row  is  unchanged  since  node  2  is  not  linked  to  any 

crack  node.  If  the  beams  had  different  properties,  then  each  term  in  U 
would  be  multiplied  by  the  corresponding  value  except  the  terms 
Ma  *  I/884.  In  the  terms,  the  sensitivity  numbers,  e^,  are  sub¬ 
scripted  by  the  original  number  of  the  column  or  girder  on  which  the 
corresponding  cracks  are  located.  Hence,  04,  for  instance,  refers  to  the 
crack  located  on  the  column  numbered  4  in  the  uncracked  frame  (Fig.  18). 
It  is  also  noted  that  Q4,  W4,  P4,  etc.  for  the  cracked  frame  have  differ¬ 
ent  values  from  the  ones  in  the  no-crack  case.  The  Z  and  X  matrices  of 

m  m 

the  cracked  frame  are  subsequently  conjugated  as: 


where  in  Equation  (78)  and  Equation  (81)  are  basically  the  same 

since  there  is  no  crack  on  the  first  wall.  The  terms  -Wu'  and  Wu'  in 
(2) 

X  In  Equation  (78)  are  replaced  by  zeros  in  Equation  (81),  since  nodes 
3  and  4  are  no  longer  linked  directly.  As  the  column  under  the  node  4  is 
cracked,  the  node  4  is  now  linked  to  the  node  7'  by  the  column  segment 

4.  Therefore,  appears  on  the  second  row  and  second  column  (which 
corresponds  to  y7)  of  F  .  The  third  row  of  X  reflects  the  continuity  of 

shear  at  the  column  crack  in  Figure  18b.  The  lower  crack-node  7  is 
linked  to  the  node  3  by  the  column  segment  7  and  the  upper  crack-node  7' 
is  linked  to  the  node  4  by  the  segment  4.  Q'  and  W*  values  are  accord¬ 
ingly  placed  on  the  third  row.  The  construction  of  the  matrices  for  the 
cracked  frame  is  now  complete. 


Figure  19  illustrates  the  relative  changes  in  the  first  three  natu¬ 
ral  frequencies  of  a  two-story  single-span  frame  with  a  single  crack  of 
8  *  .174.  The  frame  has  uniform  properties  (i.e.,  n^  »1  for  all  1,j) 
and  all  six  of  the  beam  elements  are  of  equal  length,  (I.e.,  s-f  are  the 
same  for  all  the  six  elements).  Due  to  symmetry,  a  crack  on  element  3  or 


4  will  yield  the  same  R  curves  as  a  similar  crack  on  element  1  or  2, 
respectively  (Fig.  18a).  The  first  three  characteristic  values  for  such 
a  frame  with  no  crack  are  8^  *  1.2240,  8^a  2.2177,  a  3.2885. 
Similar  to  the  guidelines  established  previously  for  damage  diagnosis  on 
a  cantilever  beam,  the  following  conclusions  can  be  deduced  from  Figure 
19  for  the  specific  frame  Involved.  (Due  to  symmetry,  only  the  left  half 
of  the  frame  need  be  considered). 

a.  In  cases  when  Rj  Is  relatively  large  and  R2  and  R3  have  near¬ 
zero  values,  the  (major)  crack  Is  either  on  the  column  element  1  at  a 
distance  of  0.35-0.40  from  the  ground  or  on  the  glrddr  element  5  at  a 
distance  of  0.0-0.15  from  the  corner,  or  there  are  two  major  cracks,  one 
each  at  the  given  locations. 

b.  In  cases  when  R2  is  relatively  large  and  Rx  and  R3  are  near 
zero,  the  crack  is  on  the  column  element  2  at  0.15-0.20  distance  from 
the  lower  joint. 

c.  For  the  case  when  R3  is  relatively  large  and  Rx  and  R3  are 
small,  the  crack  is  on  element  2  at  e7  =  0.45-0.60,  or  near  the  mid-point 
of  either  of  the  girders,  or  there  may  be  a  major  crack  at  each  of  these 
locations. 

d.  When  Rx  Is  very  small  and  both  R2  and  R3  are  relatively  large 
and  have  comparable  values  (say,  different  by  no  more  than  0.01  at  the 
range  of  values  In  Figure  19),  the  crack(s)  is  (are)  on  element  1  at 
e7  *  0.55-0.65  and/or  an  element  2  at  e7  3  0.25-0.40. 

e.  In  the  case  when  R2  Is  very  small  and  R1  and  R3  are  both  large 
and  comparable  in  value,  the  crack(s)  may  be  on  element  1  at  e7  3  0.40- 
0.50  and/or  on  element  5  at  e7  3  0.25-0.35.  However,  there  are  other 
possibilities,  an  example  of  which  Is  the  presence  of  three  major  cracks, 
one  each  on  element  1  at  e7  3  0.4-0. 5,  on  element  2  at  e7  3  0.5-0. 6,  and 
on  element  5  at  e7  3  0.0-0.20. 


f.  If  all  Rj  are  comparable  In  value,  many  possibilities  exist. 

The  above  guidelines  (for  equal-legged  frame)  are  reasonable  estimates. 
They  Illustrate  the  use  of  relative  Rj  values  when  a  numerical  solution 
cannot  be  found  for  a  single  crack.  They  also  serve  the  direction  of 
probabilistic  measure  to  be  developed  later. 

Example  2. 

The  first  three  characteristic  values  6u  for  an  undamaged  two- 
story  single-span  uniform  frame  of  equal-length  elements  are  1.2240, 
2.2177,  and  3.2885.  The  following  relative  frequency  changes  have  been 
computed  from  the  measured  frequencies  of  a  damaged  frame  of  such  charac¬ 
teristics:  Rxs  0.011,  R2  s  0.025,  R3  *  0.005.  The  damage  will  be 
estimated  based  on  these  relative  changes. 

The  characteristic  values  for  the  damaged  frame  are  computed  from 
Equation  (54).  The  first  two  of  these  values  are  used  In  the  Inverse 
problem  with  k=l  (section  3.2);  that  is,  a  solution  for  a  single  (major) 
crack  Is  sought.  Four  possible  damage  pairs  are  thus  computed;  namely, 

(e7,  ax)  =  (0.811,  0.278) 

(e7,  e6)  *  (0.094,  0.191) 

and  due  to  symmetry 

(e7,  03)  «  (0.811,  0.278) 

(e7,  0g)  *  (0.906,  0.191) 

where  e7  *  L7/L  Indicates  relative  crack  location  measured  r-om  the  lower 
end  of  each  column  or  the  left-hand  end  of  each  girder.  The  subscripts 
Neq"  for  equivalent  which  were  used  in  Chapter  3  have  been  dropped  here. 
The  subscript  on  0  indicates  the  structure  element  on  which  the  crack  Is 
located.  The  major  crack  is  thus  either  on  one  of  the  first  story 
columns  (l.e.,  on  column  1  or  3)  or  on  the  second  floor  girder  (I.e., 
beam  number  6)  (Fig.  18a).  To  decide  where  the  actual  crack  Is  located, 
(e7,  @x )  and  (e7,  @6)  are  each  used  separately  In  the  forward  problem  to 
compute  the  relative  reduction  In  the  third  natural  frequency  which  would 
be  caused  by  each  crack  alone.  The  computed  values  are  0.017  and  0.002, 


respectively.  The  second  value  is  closer  to  the  actual  measured  R3  value 
of  0.005.  Thus,  it  may  be  concluded  that  the  crack  is  on  the  second 
floor  girder,  i.e.,  on  element  6.  The  fact  that  the  actual  R3  value  is 
greater  than  the  computed  value  may  be  due  to  a  minor  crack  on  the  frame, 
for  example,  near  the  middle  of  girder  5  or  6  (Fig.  19).  It  must  be 
stressed,  however,  that  one  major  crack  is  assumed  to  exist  in  the  above 
computations.  The  same  Rj  values  may  be  produced  by  a  combination  of 
minor  cracks  or,  for  instance,  by  two  major  cracks,  one  on  column  1  at 
0.40-0.50,  the  other  on  girder  6  at  0.05-0.10  (Fig.  19).  In  the  latter 
case,  the  crack  on  column  1  would,  by  itself,  produce  almost  no  change  in 
the  second  natural  frequency  while  rendering  small  changes  (<0.003)  in 
the  fundamental  and  third  frequencies.  The  crack  on  girder  6  would,  on 
the  other  hand,  produce  negligible  change  in  the  third  frequency  while 
causing  changes  greater  than  12  in  the  other  two  frequencies.  The 
effects  of  cracks  are,  by  no  means,  linearly  additive.  Qualitatively, 
however,  a  certain  crack  contributes  most  to  the  decrease  of  that  fre¬ 
quency  which  would  experience  the  largest  decrease  if  that  particular 
crack  were  the  only  one  present  on  the  structure.  This  trend  was  Illus¬ 
trated  by  Example  5  on  a  cantilever  beam  in  the  previous  chapter.  The 
uncertainty  of  damage  diagnosis  emerges  from  the  measurability  of  only  a 
few  frequencies  of  a  structure. 

Damage  diagnosis  can  be  highly  sensitive  to  Rj  values  at  certain 
locations  of  the  Rj  curves.  If,  for  example,  Rj  is  measured  to  be 
0.010  Instead  of  0.011  In  this  example,  then  the  possible  damage  pairs 
will  be 

(e7,  81)  *  (0.795,  0.299)  ;  R3eq  =  0.018 

(e7,  06 )  *  (0.295,  0.963)  ;  R3eq  *  0.094 

together  with  their  symmetrical  counterparts. 

There  is  a  considerable  change  In  the  second  pair  while  the  first  Is 
affected  only  slightly.  In  this  case,  the  first  pair  would  be  selected 

as  the  likely  diagnosis  since  R3@q  is  closer  to  the  measured  value  of 


5.  CONCLUSION 

In  a  previous  report,  damage  diagnosis  in  simple  structures  was 
studied  [1].  In  this  report,  analogy  between  beams  and  simple  electrical 
circuits  has  been  utilized  to  extend  the  modal  method  of  diagnosis  to 
multiple  cracks  and  more  complex  frame  structures.  Computation  of  modal 
frequencies  of  an  n-story  m-span  frame  with  k  cracks  (the  forward  prob¬ 
lem)  requires  inverting  a  matrix  (or,  equivalently,  solving  a  linear 
system)  of  order  (n+k)  and  finding  zeros  of  a  deteminant  of  order 
[n(m+l)  +  2k].  In  finite  element  methods,  using  generally  20  elements 
per  basic  beam  segment,  there  are  slightly  less  than  20n(2m+l)  nodes  each 
with  two  degrees  of  freedom.  The  size  of  computation  savings  with  analog 
circuits  is  almost  40-fold.  However,  using  circuit  analogy  in  determin¬ 
ing  the  location  and  intensity  of  a  crack  (the  inverse  problem  of  diagno¬ 
ses  with  k=l)  may  involve  n(2m+l)  computations  of  zeros  of  a  determinant 
of  order  [n(m+l)  +  2]  while  inverting  a  matrix  of  order  (n+1)  each  time, 
since  the  crack  can  be  on  any  one  of  the  n(2m+l)  elements  of  the  frame. 

On  the  other  hand,  modal  frequencies  of  a  cantilever  beam  with  k  cracks 
can  be  computed  with  analog  T  circuits  by  inverting  a  matrix  and  finding 
zeros  of  a  determinant,  both  of  the  order  (k+1). 

It  is  shown,  in  the  report,  that  closely  spaced  multiple  cracks  can 
be  effectively  represented  by  a  single  crack  for  which  the  sensitivity 
number  8  is  approximately  equal  to  the  sum  of  the  individual  sensitivity 
numbers.  Uncertainty  exists,  however,  in  diagnosing  such  a  damage  as  to 
whether  there  is  only  one  major  crack  or  several  closely  spaced  minor 
cracks.  Nevertheless,  location  of  the  damage  can  be  identified  quite 
accurately. 

In  the  case  of  a  structure  with  several  cracks,  only  one  of  which  is 
severe,  it  is  possible  to  diagnose  the  major  crack.  The  contributions  of 
minor  cracks  to  frequency  decreases  is  small  compared  to  that  of  the 
severe  crack.  The  relative  frequency  changes,  Rj,  in  such  a  situation, 
exhibit  a  trend  similar  to  that  which  would  be  observed  were  the  major 
crack  the  only  crack  present  on  the  structure.  If  there  Is  more  than  one 
major  crack,  not  closely  spaced,  on  the  structure,  it  is,  in  general,  not 
possible  to  diagnose  deterministically  the  damage  with  only  three  Rj 
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values  known,  although  the  minimum  number  of  cracks  may  be  estimated. 
Currently,  the  effect  of  damping  is  being  investigated  as  a  means  to 
render  diagnosis  of  such  damage.  Curves  of  relative  frequency  changes, 

Rj  versus  crack  location  yield  valuable  Information  regarding  damage 
diagnosis.  It  may  be  psssible  to  estimate  damage  location,  provided 
there  is  only  one  major  crack,  by  only  looking  at  the  relative  magnitudes 
of  frequency  decreases  and  using  the  relative-change  curves.  The  process 
may  be  tedious  for  more  complex  structures.  Nevertheless,  such  curves 
can  be  of  qualitative  help  to  the  numerical  results.  However,  when 
cross-over  occurs  in  the  modal  frquencies,  Rj  curves  cannot  be  used 
unless  the  correspondence  between  the  pre-and  post-damage  frequencies  Is 
known.  The  approach  described  in  Section  3.2  eliminates  the  effect  of 
the  cross-over  problem  on  diagnosis  of  a  major  crack  by  eliminating  the 
need  to  know  this  correspondence. 

In  diagnosing  damage  in  a  structure,  given  three  post-damage  charac¬ 
teristic  values,  e,  the  solution  is  sought  for  a  single  crack  since  there 
are  not  enough  equations  to  solve  for  a  larger  number  of  cracks.  It  is 
shown  that  the  same  a  values  can  be  produced  by  one  crack  or  a  combina¬ 
tion  of  several  cracks.  Consideration  must  then  be  given  to  whether  It 
is  more  likely,  during  a  strong  motion  event,  to  produce  a  few  major 
cracks  or  several  minor  cracks.  In  solving  for  a  single  crack,  (i.e., 
for  eeq  and  9eq)  two  measured  a  values  are  actually  used  first.  The 
third  measured  a  value  is  then  compared  with  the  computed  counterpart, 
which  would  be  produced  by  the  crack  with  the  characteristic  pair 
(eeq,aeq)*  It  is  shown  that  the  two  a  values  may  not  match  closely 
even  though  the  diagnosis  is  accurate  {that  is,  e  and  9gq  identify  the 
major  crack  correctly),  the  discrepancy  possibly  being  due  to  the  pre¬ 
sence  of  other  minor  cracks.  If  the  discrepancy  is  large,  then  the  solu¬ 
tion  is  rejected  with  the  conclusion  that  there  is  more  than  one  major 
crack.  On  the  other  hand,  the  choice  of  the  pair  of  measured  $  values 
to  be  used  in  computing  e  and  9ec)  may  affect  the  diagnosis.  In  some 
cases,  choosing  the  two  characteristic  values  which  correspond  to  the 
largest  two  Rj  values  leads  to  the  correct  diagnosis.  It  is  not  clear, 
however,  what  the  right  choice  is  in  each  case. 


Finally,  for  frame  structures,  the  change  In  frequencies  due  to 
cracks  Is  relatively  small.  Furthermore,  damage  diagnosis  can  be  highly 
sensitive  to  Rj  values  In  certain  cases.  Therefore,  accuracy  In  field 
measurement  of  the  structural  response  Is  of  extreme  Importance. 


APPENDIX 


1.  T-CIrcult  Equations: 


£i  Z  Z  E, 


+ 1  - 


i  + 


Ci  Q  ;; 


In  this  part,  the  circuit  equations  which  relate  the  Input-output 
voltages  (y x 1 ,  y2* )  of  a  T  circuit  to  the  input-output  currents  (M1,  M2) 
are  developed.  The  terras  Input  and  output  refer  to  the  left  and  right 
ports  of  the  circuit  In  the  figure  and  are  used  In  a  relative  sense.  In 
Figure  2  in  the  text,  the  currents  and  M2  were  shown  at  the  branches 
of  the  circuit  (as  input  and  output  currents)  whereas  in  the  above  figure 
they  are  shown  as  mesh  currents.  In  principle,  there  Is  no  difference 
between  the  two.  However,  the  latter  Is  conceptually  more  helpful  In 
developing  the  circuit  equations. 

The  two  sides  of  the  circuit  in  the  figure  can  be  considered  as 
meshes.  The  sum  of  the  voltages  around  each  mesh  Is  equal  to  zero 
(Kirchhoff's  voltage  law).  Hence, 

Ex  +  (Z  +  Z0 )M i  -  Z0M2  -  yj'  »  0 
-  E2  +  (Z  +  Z0)M2  -  ZqM!  +  y2‘  *  0 

where  E  and  Z  are  the  voltage  source  and  impedance,  respectively. 
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These  equations  can  be  solved  for  the  voltages;  that  is, 

Yi'  *  (Z  +  Z0)Mj  -  Z0M2  +  Ej 
y2'  *  Z0MX  -  (Z  +  Z0)M2  +  E2 

which  are  the  same  as  Equations  (8)  in  the  text. 

2.  n-Circult  Equations: 
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In  this  circuit,  there  are  two  nodes,  with  node  voltages  and 
y2',  and  a  reference  node,  y0‘  which  Is  assigned  a  zero  voltage.  G  and 
G0  are  conductances  and  ^  and  I2  are  current  sources.  The  sum  of  the 
currents  entering  each  node  Is  zero  (Kirchhoff's  current  law).  Hence, 
for  nodes  1  and  2 

Mx  -  Ij  -  G(y1‘  -  y0* )  -  G„(y -  y2')  »  0 
-  M2  +  I2  -  G(y2*  -  y0* )  -  G0(y2'  -  y^)  *  0  . 

The  term  Goty^  -  y2‘),  for  Instance,  denotes  the  current  flowing  through 
the  conductance  G0  from  the  node  1  to  the  node  2.  There  is  a  minus  sign 
In  front  of  the  term  since.  In  the  first  equation,  currents  entering  the 
first  node  are  taken  to  be  positive.  The  equations  can  be  solved  for  the 
currents  Mj  and  M2;  hence,  with  y0'  *  0, 

Mi  *  (G  +  Gq )y j *  -  G(jy2'  +  Ii 
M2  *  Ggy i'  *  (G  +  G o )y 2*  +  ^2 

which  are  the  same  as  Equations  (21)  In  the  text. 


9 

\ 


r 


67 


